1 991AJ . . . . 102 . 151 OH 


THE  ASTRONOMICAL  JOURNAL  VOLUME  102,  NUMBER  4  OCTOBER  1991 

THE  MOTION  OF  MARS’  POLE.  I.  RIGID  BODY  PRECESSION  AND  NUTATION 

James  L.  Hilton 

U.S.  Naval  Observatory,  34th  Street  and  Massachusetts  Avenue  NW,  Washington,  DC  20392 
Received  5  March  1991;  revised  20  June  1991 

ABSTRACT 

The  total  precession  in  longitude  for  a  solid-core,  rigid  Mars  model  is  found  to  be 
AT  =  —  7"296  +  0"021  yr“ 1 .  This  precession  includes  the  motion  resulting  from  the  extremely  long- 
period  ( >  10  000  yr)  nutation  components  that  result  from  changes  in  Mars’  orbit  with  time.  These 
same  very  long-period  nutation  components  also  contribute  a  “precession”  in  latitude  of 
Ac  =  0"4255  +  0"0012  yr  - 1 .  The  nutation  components  in  longitude  with  amplitudes  larger  than  0"001 
consist  of  eight  solar  nutation  terms  and  one  nutation  each  contributed  by  the  motions  of  the  nodes  of 
Phobos  and  Deimos  and  a  single  nutation  from  the  direct  gravitational  torque  of  Jupiter.  The  nutation 
components  in  obliquity  consist  of  six  solar  terms  and  one  nutation  term  each  from  the  motion  of  the 
nodes  of  Phobos  and  Deimos.  The  amplitudes  and  periods  of  the  solar  nutation  terms  are  in  agreement 
with  the  nutation  found  by  Reasenberg  &  King  [J.  Geophys.  Res.,  84,  6231  (1979)].  The  solar  preces¬ 
sion  rate  also  agrees  with  Reasenberg  and  King’s  value  with  the  two  small  corrections  resulting  from 
the  change  in  the  orientation  of  Mars’  orbit  with  time  added  to  it.  A  single  nutation  in  longitude  driven 
by  Jupiter  is  the  only  significant  planetary  contribution  at  the  milliarcsecond  level. 


1.  INTRODUCTION 

Differences  between  the  observed  forced  nutation  compo¬ 
nent  amplitudes  of  the  Earth  and  predicted  nutation  ampli¬ 
tudes  for  a  rigid  Earth  are  a  result  of  differences  between  the 
theoretical  rigid  structure  of  the  Earth  used  in  older  models 
and  the  actual  elastic  Earth  with  a  liquid  core.  However, 
except  for  the  period  of  the  Chandler  wobble,  the  observa¬ 
tions  of  the  motion  of  the  Earth’s  pole  were  not  accurate 
enough  to  observe  the  effects  of  the  elastic,  liquid  core  Earth 
until  the  last  30  years.  In  more  recent  works,  such  as  those  by 
Wahr  ( 1981a, b),  the  nutation  resulting  from  an  elastic,  liq¬ 
uid  core  Earth  are  modeled  as  perturbations  of  the  rigid 
Earth  model  nutation.  This  is  the  approach  adopted  for  the 
1980 IAU  Theory  of  Nutation  to  determine  the  amplitude  of 
the  various  nutational  elements  (Kaplan  1981).  These  per¬ 
turbations  result  in  modifications  to  the  nutation  amplitudes 
for  the  Earth  from  about  1%  to  0.01%  of  the  theoretical 
rigid  nutation  amplitudes  or  about  0"0 1 9  for  the  largest  term 
in  the  series.  Since  the  1950s  the  improvement  in  the  mea¬ 
surement  of  the  motion  of  the  Earth’s  pole  have  made  its 
observation  a  powerful  probe  of  the  structure  of  the  Earth. 
The  recent  increase  in  the  quality  of  the  data  available  for  the 
precession  and  nutation  for  the  Earth  has  started  a  reevalua¬ 
tion  of  nutation  theory  for  both  rigid  and  elastic  Earth  mod¬ 
els  to  improve  their  accuracies  (e.g.,  Kinoshita  &  Souchay 
1990;  Zhu  etal.  1990).  The  most  recent  work,  such  as  done 
with  VLBI,  has  made  the  information  on  the  Earth’s  interior 
obtained  from  observations  of  precession,  nutation,  and  po¬ 
lar  motion  of  the  same  quality  as  can  be  determined  from 
seismometry  (Melchiorl986).Thismethodcanbe  used  as  a 
probe  of  the  structure  of  other  planets  in  the  solar  system 
provided  that  it  is  possible  to  observe  the  planet’s  orientation 
in  space  with  high  enough  accuracy.  This  requirement  can 
definitely  be  fulfilled  for  only  one  other  planet  in  the  solar 
system,  Mars. 

Mars,  the  fourth  planet  from  the  Sun,  is  the  third  in  size  of 
the  terrestrial  planets  of  the  solar  system.  Its  orbit  brings 
Mars  to  within  0.53  AU  of  the  Earth.  Its  surface  is  not  cov¬ 
ered  by  clouds  like  Venus  and  the  gaseous  giant  planets,  and, 
unlike  Mercury,  it  is  observable  away  from  the  Sun’s  glare. 
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These  properties  make  Mars  the  easiest  planet  in  the  solar 
system,  aside  from  the  Earth,  to  observe.  Despite  its  small 
size  [equatorial  radius  3393.4  km  and  mass  6.42 X 1023  kg; 
Astronomical  Almanac  for  the  Year  1990  (1989)],  Mars 
shows  many  Earthlike  properties  in  its  physical  ephemeris. 
Mars  has  a  sidereal  day  of  1 .0260  Earth  solar  days,  the  incli¬ 
nation  of  its  equator  to  the  plane  of  its  orbit  of  25°20,  and  a 
geometric  oblateness  of  0.0052  in  comparison  to  the  Earth’s 
geometric  oblateness  of  0.0034. 

It  is,  however,  a  mistake  to  assume  that  Mars  has  an  inter¬ 
nal  composition  that  is  the  same  as  the  Earth.  First,  Mars  has 
a  mean  density  that  is  only  61.9%  of  the  mean  density  of  the 
other  terrestrial  planets. 

The  inertia  ratio  is  a  measure  of  a  planet’s  central  conden¬ 
sation  defined  by 

q=C/^R2,  (1) 

where  C  is  the  principal  moment  of  inertia  about  the  polar 
axis,  is  the  mass  of  the  planet,  and  R  is  the  equatorial 
radius  of  the  planet.  The  estimates  of  q  for  Mars  range  from 
0.3654  (Reasenberg  1977)  to  0.3452  (Bills  1989)  in  com¬ 
parison  to  0.3335  for  the  Earth  and  0.4  for  a  homogeneous 
sphere.  The  large  value  for  the  inertia  ratio  means  that  Mars 
is  less  centrally  condensed  than  the  Earth,  but  it  is  also  small 
enough  to  show  that  Mars  is  not  completely  homogeneous. 
Mars  may  have  a  very  small  external  magnetic  field,  the 
strength  of  which  has  not  yet  been  accurately  measured,  but 
it  has  been  established  that  Mars’  magnetic  field  is  much 
smaller  than  the  Earth’s.  Other  information  on  the  interior 
structure  of  Mars  is  very  sketchy.  The  only  structural  infor¬ 
mation  that  is  certain  is  Mars’  radius  and  its  mean  density.  In 
addition  the  structure  of  the  surface  gravity  field  to  twelfth 
degree  and  order  has  been  published  by  Christensen  &  Bal- 
mino  ( 1979),  but  the  higher  order  structure  of  the  field  is 
only  preliminary. 

The  seismic  experiments  included  in  the  Viking  probes 
were  designed  only  as  a  preliminary  survey  to  determine 
what  instruments  would  be  necessary  on  future  missions  to 
Mars.  Also,  the  seismic  package  functioned  properly  on  only 
one  of  the  two  landers  (Anderson  et  al.  1977). 
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The  less  extreme  central  condensation  of  Mars  and  its 
small  magnetic  field  are  the  reasons  that  some  of  the  existing 
models  of  Mars  picture  a  planet  with  a  solid  core  (Binder  & 
Davis  1973 )  rather  than  a  liquid  core  as  the  Earth  has.  Most 
models,  however,  do  use  liquid  cores,  but  the  composition  of 
the  core  is  very  much  in  question.  Johnston  &  Toksoz 
( 1977)  and  Okal  &  Anderson  ( 1978)  discuss  the  merits  of 
various  possible  liquid  core  configurations.  The  two  main 
types  of  core  used  in  the  models  are  either  small  and  formed 
from  molten  Fe  or  large  and  containing  a  molten  iron  com¬ 
pound  such  as  FeS.  If  the  core  is  composed  of  Fe,  it  needs  to 
be,  proportionately,  much  smaller  than  the  Earth’s  core  to 
be  consistent  with  the  smaller  central  condensation  of  Mars 
and  to  account  for  the  small  magnetic  field.  A  core  com¬ 
posed  of  FeS  would  be  larger  than  a  Fe  core  since  FeS  is  less 
dense.  Also,  FeS  is  not  as  good  an  electrical  conductor  as  Fe, 
so  the  small  magnetic  field  is  not  a  consideration.  The  range 
of  models  investigated  by  Okal  &  Anderson  ( 1978 )  included 
cores  with  radii  from  20%  to  60%  of  Mars’  radius  and  con¬ 
tained  7%  to  31%  of  its  total  mass.  The  key  to  determining 
the  core’s  composition  is  determining  its  mean  radius  which, 
along  with  the  planetary  inertia  ratio,  will  determine  its 
mean  density. 

By  using  the  Earth  as  an  example  for  how  the  motion  of 
the  pole  is  affected  by  a  planet’s  structure  and  keeping  in 
mind  the  known  and  possible  differences  between  Mars  and 
the  Earth,  it  is  possible  to  use  the  motion  of  Mars’  pole  as  a 
probe  of  its  interior  structure.  The  information  needed  to 
discriminate  between  the  various  models  of  the  Martian  inte¬ 
rior  can  be  obtained  from  a  few,  highly  accurate  planetary 
orientation  sensors  as  opposed  to  a  large  net  of  seismometers 
needed  for  a  classic  determination  of  the  planetary  interior. 

Mars  has  been  the  subject  of  most  of  the  attempts  to  deter¬ 
mine  the  precession  and  nutation  of  other  planets.  One  of  the 
first  attempts  to  calculate  Mars’  precession  was  by  Struve 
(1898)  who  obtained  a  value  of  —  7"07  per  terrestrial  year. 
More  recent  attempts  to  determine  the  precession  of  Mars 
theoretically  are  those  of  de  Vaucouleurs  (1964)  with  a 
precession  rate  of  —  7"07  per  year  and  Reasenberg  &  King 
(1979)who  determined  a  precession  rate  of  — 7"575  per 
year  along  with  nine  nutation  components  in  longitude  and 
seven  nutation  components  in  obliquity.  These  nutation 
components  are  the  result  of  the  direct  solar  torque  and  long- 
period  changes  in  Mars’  orbit  only.  The  nutation  compo¬ 
nents  range  in  amplitude  from  1  "097  to  0"0004  with  periods 
ranging  from  1 14.5  days  to  686.9  days.  Reasenberg  &  King 
also  made  a  first  study  of  the  effects  of  other  gravitational 
sources  on  the  above  terms.  However,  some  sources  for  mo¬ 
tion  of  the  pole  have  not  been  taken  into  account,  such  as  ( 1 ) 
tidal  distortion,  that  is,  planetary  elasticity,  (2)  diurnal  po¬ 
lar  motion  due  to  tidal  distortion,  ( 3 )  regional  tectonocism 
which  may  be  responsible  for  the  buildup  of  the  Tharsis  re¬ 
gion  of  Mars  (Reasenberg  1977)  and  the  formation  of  the 
Valles  Marineris  (Blasiusef  a/.  1977),  and  (4)  the  torques  of 
the  Martian  satellites  on  the  precession  and  nutation  of 
Mars.  Although  both  satellites  are  small  and  in  nearly  equa¬ 
torial  orbits,  they  are  close  enough  to  Mars  that  preliminary 
calculations  indicate  that  Phobos  has  the  potential  for  pro¬ 
ducing  a  precession  of  at  least  a  few  tenths  of  an  arcsecond 
per  year.  Ward  (1974)  and  Borderies  (1980)  have  done 
work  on  the  long  period  nutation  of  Mars. 

A  preliminary  study  of  the  motion  of  Mars’  pole  from 
Viking  data  has  been  carried  out  by  Borderies  etal.  ( 1980) . 
This  paper  uses  the  data  from  the  Martian  ranging  and 


Doppler  shift  information  from  the  first  nine  months  of  the 
Viking  landers  on  Mars.  The  results  from  these  preliminary 
measurements  show  ( 1 )  there  are  definite  systematic  residu¬ 
als  left  from  their  solution  for  the  planetary  motion  and  (2) 
the  data  do  not  provide  a  long  enough  baseline  for  a  defini¬ 
tive  measurement  of  the  motion  of  Mars’  pole. 

The  object  of  this  paper  is  to  determine  the  theoretical 
basis  for  the  rigid-body  precession  and  nutation  of  Mars  in¬ 
cluding  the  contributions  of  Phobos,  Deimos,  and  the  plan¬ 
ets.  A  future  paper  will  take  the  rigid-body  framework  devel¬ 
oped  here,  and  apply  it  to  an  elastic  planet.  Models  of  Mars 
using  either  solid  or  liquid  core  configurations  will  show 
how  the  radius  of  the  core  and  the  core’s  state  affect  the 
precession,  nutation,  and  other  elements  of  the  motion  of 
Mars’  pole.  This  pair  of  papers  will  then  provide  a  basis  for 
the  first-order  determination  of  Mars’  structure  based  on  the 
change  in  its  spatial  orientation  with  time. 

2.  THE  EQUATION  OF  MOTION 

The  gravitational  potential  of  a  test  particle  at  r  by  an 
extended  mass ..//  is  expressed  by 

<D(r)  Jp,(rT')fi^ 

+  ±  jp2(r-r')dm  +  -■),  (2) 

where  r'  is  the  position  of  an  infinitesimal  mass  element  dm 
within  the  extended  mass  and  P„  (r*r')  are  the  Legendre 
polynomials  of  order  n.  The  geometric  algebra  form  of  the 
Legendre  polynomial  uses  the  dot  product  of  two  vectors 
rather  than  the  usual  cosine  of  the  angle  between  the  vectors 
for  the  free  parameter. 

From  Appendix  A  the  first-order  equation  for  the  torque 
of  a  point  mass  m  on  a  spheroidal  planet  is 


where  m  is  the  mass  of  the  perturbing  body,  r  is  the  distance 
from  the  center  of  mass  of  the  planet  to  the  perturbing  body, 
r  is  the  unit  vector  pointing  toward  the  perturbing  body  from 
the  center  of  mass  of  JS ,  e  is  the  unit  vector  along  the  axis  of 
the  greatest  principal  moment  of  inertia,  J2  is  the  dynamical 
form  factor,  and  R  is  the  equatorial  radius  of  the  planet. 

The  equation  for  the  forced  motion  of  the  pole  of  a  spher¬ 
oidal  body  is  found  to  be 

(  —  A  sin  e,  —  (3, 0) 

=  ^.(sin  8  cos  8  sin  a, sin  8  cos  8  cos  a, 0),  (4) 

r*n  q 

where  A  is  the  time  derivative  of  the  motion  in  longitude  with 
respect  to  a  given  inertial  coordinate  plane,  j3  is  time  deriva¬ 
tive  of  the  latitude,  a  is  the  longitude  of  the  perturbing  body 
with  respect  to  the  equator  of  the  planet,  and  8  is  the  latitude 
of  the  perturbing  body  with  respect  to  the  planet’s  equator. 

3.  SOLAR  PRECESSION 

For  the  solar  precession  and  nutation  two  additional  sim¬ 
plifications  can  be  made.  First,  the  inclination  of  the  orbital 
plane  of  the  Sun  about  Mars  and  the  orbital  plane  of  Mars 
about  the  Sun  are  identical,  so  e  =  e0.  Second,  the  mean  orbi¬ 
tal  plane  at  the  chosen  epoch  is  identical  to  the  reference 
plane,  that  is  /?  —  0,  so  cos  P  =  1  and  sin  /3  =  0. 
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3. 1  Precession  and  Nutation  for  a  Circular  Orbit 

In  Eq.  (4)  the  angle  e  varies  only  a  small  amount  over 
relatively  short  periods  of  time  (~104  yr).  Ward  (1974) 
and  Borderies  ( 1980)  showed  that  the  obliquity  of  Mars 
oscillates  about  4°  on  a  time  scale  of  1 . 2  X 1 0s  yr  and  about  7° 
over  1.2  X  106  yr.  These  changes  in  inclination  are  mainly  a 
result  of  changes  in  the  orbital  plane  of  Mars  rather  than 
changes  in  the  orientation  of  the  rotation  axis.  Thus,  for  the 
present  purposes,  e  is  considered  constant.  The  effect  of  the 
change  in  the  obliquity  of  Mars  on  its  precession  and  nuta¬ 
tion  will  be  covered  in  Sec.  3.3  and  Appendix  D.  Therefore, 
for  a  circular  orbit,  the  only  quantity  that  changes  with  time 
on  the  right-hand  side  of  Eq.  (4)  is  the  longitude  A.  The 
zeroth-order  approximation  for  A  is  a  monotonically  in¬ 
creasing  function  of  time: 

A(t)  -nt  +  A0,  (5) 

where  A0  is  the  ecliptic  longitude  at  some  initial  epoch  t0. 
Using  Eq.  (5)  in  Eq.  (4)  gives 

A  sin  e  = - —  —  ^  sin  e  cos  e  sin 2(nt  +  A0)  (6) 

Cl  q  r 

and 

=  —  ^-sin(nt  +  A0)sinecos(nt  +  A0)  .  (7) 

Cl  q  r 

These  two  equations  for  A  and  (3  are  integrated  to  determine 
both  the  precession  in  longitude  and  the  first  nutation  terms 
in  longitude  and  latitude  that  arise  from  a  circular  orbit  of 
Mars  about  the  Sun.  Equation  (6)  is  divided  by  sin  e  and 
integrated  over  the  orbital  period.  The  change  in  longitude 
of  the  pole  as  a  function  of  time  is 

A  =  — 2^-  CL  cos  J_ L  t  -  -2-  sin  V,nt  +  A0))  +  C. 

Cl  q  r  \  2  4  n  J 

(8) 

Integrating  Eq.  (27)  gives 

jg=  ~~  ^  —  ^-sin  e— j—  cos  2(«r  +  A0)  +  C'.  (9) 

Cl  q  r  4  n 

Using  the  definition  of  the  Euler  angles  Ae  =  —  A/?  and 
A ip  =  —  A A.  Thus,  the  negatives  of  Eqs.  (8)  and  (9)  give 
the  precession  in  if)  and  the  circular  orbit  nutation  terms  in 
both  e  and  ip. 

For  the  solar  precession  and  nutation  of  Mars,  the  values 
for  the  quantities  needed  to  evaluate  Eqs.  (8)  and  (9)  are 

G=  (6.672  +  0.001  )X  10"  11  m’kg-'s”2 

(Kaplan  1981), 

^  =  (1.9891  ±0.0001  )X1030  kg  (Kaplan  1981), 
r=  (2.279  390  77X  10n  ±  1500)m 
(Kaplan  1981;  Bretagnon  1982), 

J2  =  (1.959  0468  ±0.000  138  70)Xl0”3 
(Christensen  &  Williams  1978), 
q  =  0.3654  ±  0.0010  (Reasenberg  1979), 
fl  =  (7,088  2181Xl0-5±  lXlO-12)s-‘ 

(Mayo  etal.  1977), 

e  =  (25° 20  ±  0!02)  (Sinclair,  1972). 
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The  exact  value  for  the  inertia  ratio  q  is  still  unknown  as 
pointed  out  by  Bills  ( 1989).  The  value  of  q  is  derived  from 
models  of  Mars’  gravity  field.  These  gravity  field  models 
incorporate  an  assumed  orientation  for  the  nonhydrostatic 
component  of  Mars’  mass  distribution  with  respect  to  its 
equator.  Reasenberg’ s  (1977)  value  finds  that  the  pole  of  the 
nonhydrostatic  component  of  Mars’  structure  is  oriented 
along  the  pole  of  Mars’  greatest  principal  moment  of  inertia. 
Bills,  however,  argues  based  on  statistics  that  the  nonhy¬ 
drostatic  component  of  Mars’  mass  distribution  is  more  like¬ 
ly  to  have  the  pole  of  Mars’  nonhydrostatic  component  in  the 
midlatitude  regions.  For  this  orientation  of  the  nonhydrosta¬ 
tic  component  the  best  value  for  Mars’  inertia  ratio  is 
q  =  0.345.  However,  the  two  largest  visible  nonyhydrostatic 
components  of  Mars,  the  Valles  Marineris  and  the  Tharsis 
uplift,  are  both  centered  near  the  equator,  and  Mars,  unlike 
the  Earth,  has  a  definite  correlation  between  its  low-order 
gravitational  field  components  and  its  terrain  features  (Lor- 
ell  et  al.  1 972 ) .  Other  arguments  favoring  Reasenberg’s  iner¬ 
tia  ratio  are  presented  by  Kaula  et  al.  ( 1 9  8 9 ) .  Therefore  the 
value  for  the  inertia  ratio  used  here  is  0.3654  (Reasenberg 
1977). 

Using  the  above,  constants  for  Mars,  the  precession  of 
Mars,  the  first  term  in  Eq.  (28)  is 

ip=  -  (1.1504±0.0032)X10-12  s~‘ 

=  -  (7"488  ±  0"021 )  yr”1. 

This  value  for  the  mean  precession  of  Mars’  pole  compares 
very  favorably  with  previously  derived  values  for  the  polar 
precession.  By  far,  the  largest  source  of  uncertainty  (99.7% ) 
in  the  precession  rate  is  the  uncertainty  in  the  inertia  ratio  q. 

The  amplitude  of  the  two  circular  orbit  nutation  compo¬ 
nents  driven  by  the  torque  of  the  Sun  on  Mars  are  deter¬ 
mined  using  the  mean  orbital  motion  of  Mars  taken  from 
Bretagnon  (1982) 

n=  1.058  589  015xl0”7±  1x10” 16  rads”1. 
Substituting  this  into  Eq.  (8)  gives  a  nutation  amplitude  of 

—  SL  hCL  cos  e=  (5.434  +  0.015)  XlO”6  rad 

An  ft  q  r3 

=  17121  ±  0"003. 

Similarly,  using  the  value  for  sin  e  the  amplitude  of  the  main 
nutation  in  obliquity  determined  from  Eq.  (9): 

_zJ.^A^.sine=  (  -  2.557  ±  0.007)  X 10”6  rad 
An  Cl  q  r 

=  -  0:'527  ±  0"002. 

Both  of  these  nutation  amplitudes  are  in  good  agreement 
with  the  primary  nutation  amplitudes  determined  by  Rea¬ 
senberg  &  King  (1979)  of  1"097  in  longitude  and  0"5157  in 
obliquity. 

Now  only  the  initial  phase  A0  needs  to  be  determined.  The 
restoring  force  on  Mars  is  at  its  greatest  when  Mars  is  at  its 
solstices.  Therefore,  the  greatest  deviation  from  the  mean 
will  also  occur  at  the  solstices.  Letting  the  mean  longitude  of 
Mars  be  represented  by  Jtf  and  the  angle  between  the  perihe¬ 
lion  of  Mars  and  its  vernal  equinox  be  A,  then 

nt  +  A0  =  Jf  —  A.  (10) 

The  circular  orbit  precession  and  nutation  terms  for  Mars 
are 


©  American  Astronomical  Society 


Provided  by  the  NASA  Astrophysics  Data  System 


1 991AJ . . . . 102 . 151 OH 


1513  JAMES  L.  HILTON:  MARS’  POLE 


1513 


ip  =  —  7"488t  +  17121  cos  2(2z^  —  A), 

Se  =  0'.'527  sin  —  A).  (11) 

3.2  Components  from  the  Orbital  Eccentricity 

The  next  level  of  complication  to  add  to  the  determination 
of  the  amplitudes  of  Mars’  nutation  components  is  the 
change  in  the  torque  on  the  planet  as  a  function  of  distance 
from  the  Sun.  The  nutation  will  be  determined  by  using  a 
series  expansion  in  sine  and  cosine  functions  of  angular  vari¬ 
ables.  The  angular  variables  are  chosen  to  be  linear  functions 
of  time. 

For  the  Sun’s  torque  on  Mars,  the  angular  variables  used 
are  ff  and  A.  The  result  to  expect  from  such  an  expansion  is 
that  the  two  circular  orbit  nutation  terms  in  Eqs.  (10)  and 
(11)  above  will  be  split  into  two  series  of  nutation  compo¬ 
nents  for  the  planet.  From  Appendix  B  the  precession  and 
nutation  are 

2  =  cos  e(A0t  +  jr  — [ A„  sin  (nJf) 
a  V  n=i  nJz 

B  \ 

—  Cn  cos  (2A  +  nJ?)  ]  H - —  sin  (2A  +  nJ?) ) , 

P  -  Sin  J  £  cos  (2A  +  nJf) 

O  \n  —  o  YloZi 

+  -%-sin(2A  +  n^)),  (12) 

where  A0,  An,B„,  and  C„  are  coefficients  determined  in  Ap¬ 
pendix  B.  Inserting  the  appropriate  values  for  the  constants 
for  Mars  the  coefficient  for  the  precession  becomes 

ip  =  —  7"587  +  0"021  yr“’. 

The  coefficients  for  the  nutation  terms  larger  than  O'.'OOl 
are  given  in  Table  1.  All  of  the  nutation  amplitudes  fromC„ 
coefficients  are  smaller  than  0'.'001  so  there  are  no 
[cos  (2A  +  nJf)  ]  terms  in  longitude  or  (sin  2A  sin  nff) 
terms  in  obliquity  in  Table  1 .  As  before,  the  largest  source  of 
uncertainty  is  the  inertia  ratio  of  Mars. 


3.2.1  Changes  of  eccentricity  with  time 

Mars  is  subject  to  planetary  gravitational  torques  as  well 
as  the  torque  provided  by  the  Sun.  The  planetary  perturba¬ 
tions  cause  the  shape  and  orientation  of  Mars’  orbit  to  vary 
with  time.  The  two  changes  in  the  orbit  and  orientation 
which  are  of  concern  are  the  eccentricity  of  Mars’  orbit  and 
the  obliquity  of  Mars  to  its  orbit.  The  change  in  the  obliquity 
does  cause  small  but  significant  changes  in  the  nutation  in 
obliquity  and  will  be  treated  in  the  next  section.  In  this  sec¬ 
tion  the  effect  of  change  in  the  eccentricity  of  the  orbit  will  be 
explored. 

The  treatment  of  the  eccentricity  uses  the  Bretagnon 
(1982)  semianalytic  theory  for  Mars  and  is  divided  into  pe¬ 
riodic  and  secular  terms.  Appendix  C  shows  that  all  of  the 
periodic  components  caused  by  the  change  in  eccentricity 
are  insignificant.  Substituting  the  nonperiodic,  zeroth-order 
corrections  into  Appendix  C  the  total  precession  becomes 

ift=  -7"588  +  0':021  yr-1. 

The  change  of  0"001  yr“ 1  is  smaller  than  the  uncertainty  in 
the  precession. 

Including  the  secular  zeroth-order  part  for  eccentricity 


Table  1.  Amplitudes  for  the  rigid  body  solar  nutation  components  of 
Mars. 


Nutation  Components  in  Longitude 


Term 

Amplitude 

n 

Error 

n 

Period 

(day) 

sin  L 

-0.6343 

0.0018 

686.93 

sin  2 L 

-0.0443 

0.0001 

343.46 

sin  3 L 

-0.00405 

0.00001 

228.98 

sin  (2A  +  L) 

-0.1046 

0.0003 

686.72 

sin  (2A  +  2 L) 

1.0963 

0.0031 

343.41 

sin  (2A  +  3  L) 

0.2396 

0.0007 

228.96 

sin  (2A  +  4 L) 

0.0407 

0.0001 

171.72 

sin  (2A  +  5 L) 

0.00630 

0.00002 

137.38 

sin  (2A  +  6 L) 

0.000926 

0.000003 

114.48 

Nutation  Components  in  Latitude 

Term 

Amplitude 

Error 

Period 

c) 

c) 

(day) 

cos  (2A  +  L) 

-0.0492 

0.0001 

686.72 

cos  (2A  +  2 L) 

0.5159 

0.0014 

343.41 

cos  (2A  +  3 L) 

0.1127 

0.0003 

228.96 

cos  (2A  +  4 L) 

0.01917 

0.00005 

171.72 

cos  (2A  +  5  L) 

0.002963 

0.000008 

137.38 

L  =  Mars'  mean  longitude 

A  =  angle  between  the  perihelion  and  the  vernal  equinox  of  Mars 


from  Bretagnon’s  theory  does  cause  a  change  that  is  larger 
than  the  uncertainty  of  the  individual  terms  for  most  of  the 
nutation  terms  (Table  2 ) .  For  most  of  the  components,  how¬ 
ever,  the  change  in  amplitude  is  too  small  to  avoid  being  lost 
in  the  uncertainty  of  the  larger  nutation  components.  The 
effect  of  the  change  of  eccentricity  with  time  is  seen  as  a  small 
change  in  the  amplitudes  of  the  nutation  components. 

3.3  Changes  in  Orbital  Inclination  and  the  Ascending  Node 

A  final  source  for  solar  nutation  components  is  long  term 
change  in  the  obliquity  of  Mars’  equator  to  its  orbit.  The 
source  of  the  changes  in  obliquity  is  not  the  nutation  in  obliq¬ 
uity,  terms  which  have  periods  of  a  Martian  sidereal  year 
(686.9297  days)  or  less.  The  change  in  obliquity  is  the  result 
of  changes  in  the  orientation  of  Mars’  orbital  plane  with  time 
in  response  to  the  perturbations  of  the  other  planets  in  the 
solar  system.  Ward  ( 1974)  and  Borderies  ( 1980)  show  that 
the  variation  in  the  obliquity  of  Mars  resulting  from  changes 
in  the  orbital  plane  can  be  represented  by  a  pair  of  approxi¬ 
mately  equal  amplitude  oscillations;  the  first  oscillation  has 
a  period  of  1. 14  X 105  yr,  the  second  oscillation  has  a  period 
of  1.26X  105  yr,  and  the  combined  amplitude  is  about  10°3. 

The  equations  of  motion  for  Mars  developed  by  Ward 
( 1974)  are 

—  =  —  cos  w  sin; - 1-  sin  w —  (13) 

dt  dt  *  dt 

and 

—  (sin  i  sin  ip  cot  e  —  cos  /) 
dt  dt 


—  a  cos  e  +  cos  lAcot  e — .  (14) 

dt 

This  is  a  pair  of  coupled  first-order  equations  for  e  and  ip  as 
a  function  of  time.  These  equations  are  solved  by  making 
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Table  2.  Amplitudes  for  the  rigid  body  solar  nutation  components  of 
Mars  including  the  change  in  eccentricity  with  time. 


Nutation  Components  in  Longitude 

Term 

Amplitude 

Error 

Period 

r) 

o 

(day) 

sin  L 

-0.6357 

0.0018 

686.93 

sin  2 L 

-0.0445 

0.0001 

343.46 

sin  3 L 

-0.00408 

0.00001 

228.98 

sin  (2A  +  L) 

-0.1047 

0.0003 

686.72 

sin  (2A  +  2 L) 

1.0962 

0.0031 

343.41 

sin  (2A  +  3  L) 

0.2401 

0.0007 

228.96 

sin  (2A  +  4 L) 

0.0409 

0.0001 

171.72 

sin  (2 A  +  5  L) 

0.00634 

0.00002 

137.38 

sin  (2 A  +  6 L) 

0.000934 

0.000003 

114.48 

Nutation  Components  in  Latitude 

Term 

Amplitude 

Error 

Period 

c) 

c) 

(day) 

cos  (2A  +  L) 

-0.0493 

0.0001 

686.72 

cos  (2A  +  2 L) 

0.5158 

0.0014 

343.41 

cos  (2A  +  3  L) 

0.1130 

0.0003 

228.96 

cos  (2A  +  4 L) 

0.01925 

0.00005 

171.72 

cos  (2A  +  5  L) 

0.002982 

0.000008 

137.38 

L  =  Mars'  mean  longitude 

A  =  angle  between  the  perihelion  and  the  vernal  equinox  of  Mars 


Years  from  present 

Fig.  1.  Very  long-period  nutation  in  obliquity  over  historical  periods 
of  time.  Over  short  time  spans  ±  10  000  yr  the  change  in  the  inclina¬ 
tion  of  Mars’  pole  is  nearly  linear  with  time. 


4.  PRECESSION  AND  NUTATION  FROM  NATURAL  SATELLITES 

Mars  has  two  satellites,  Phobos  and  Deimos,  to  contribute 
to  Mars’  precession  and  nutation,  just  as  the  Moon  contrib¬ 
utes  to  the  Earth’s  precession  and  nutation.  Although  both 
of  these  satellites  are  small,  their  contributions  to  Mars’  nu¬ 
tation  may  be  fairly  large  because  their  mean  distances  from 
Mars  are  very  small  and  forced  nutation  is  a  1  /r3  effect.  The 
orbital  elements  of  the  satellites  taken  from  Sinclair  ( 1972) 
and  the  masses  taken  from  Morley  ( 1990)  are 


repeated  approximations  of  their  actual  values.  The  first- 
order  solutions  to  the  equations  are  a  very  good  approxima¬ 
tion  to  the  motion.  Ward  ( 1974)  solved  the  first-order  equa¬ 
tion  for  e  and  gives  a  first-order  estimate  to  the  first-order 
solution  for  the  motion  in  longitude  ip.  The  first-order  solu¬ 
tion  to  the  equation  of  motion  in  latitude  and  the  third-order 
estimate  to  the  first-order  solution  to  the  equation  of  motion 
in  longitude  are  developed  in  Appendix  D. 

Over  short  periods  of  time  near  the  present  era,  ±  10  000 
yr  (Fig.  1 ),  the  motion  in  obliquity  of  the  Martian  pole  as  a 
result  of  the  change  in  orientation  of  Mars’  orbit  can  be  rep¬ 
resented  by  the  linear  function 

e  =  (25°2  +  0°01 )  +  0"4255f,  (15) 

where  t  is  the  time,  in  years,  from  J2000.0.  Also,  the  zeroth- 
and  first-order  parts  of  the  first-order  solution  for  precession 
and  nutation  in  longitude  appear  as  a  contribution  to  the 
precession  of 

ipn  =  0:’2127  ±  0"0006  yr”1 
shown  in  Fig.  2. 

The  secular  portion  of  the  second-order  part  to  the  first- 
order  solution  gives  an  additional  correction  in  precession  of 

Sipn  =  0"0797  +  0:'0002  yr” '. 

Overall,  the  change  of  Mars’  orbit  with  time  causes  two 
measurable  effects  in  longitude.  ( 1 )  A  set  of  very  long-term 
nutation  components  from  the  first-order  term.  These  nuta¬ 
tion  components  are  seen  over  short  periods  of  time  as  a 
“precession”  of  0"2 137  yr  ” 1 .  (2 )  There  is  an  addition  to  the 
precession  from  the  second-order  term  of  0"0797  yr”1. 
Adding  these  terms  to  the  solar  precession  the  total  preces¬ 
sion  in  longitude  of  Mars  is 

ip  =  —  7"296  ±  0"021  yr”1. 

In  addition,  there  is  a  “precession”  in  latitude  of  0"4255 
yr”1. 


Phobos: 

m=  (1.05 +  0.1)  X1016  kg, 
r=  (9.378  +  0.001  )X  106  m, 

/  =  1°01  ±0°07, 

n  =  7196.1730  ±0.0064  rad/yr, 

A0  =  2.817  +  0.003  rad, 

N=  —  2.776  +  0.003  rad/yr, 

N—  2.65  +  0.07  rad, 

where  /  is  the  inclination  of  Phobos’  orbit  to  the  equator  of 
Mars,  n  is  the  mean  motion  of  Phobos,  A0  is  the  longitude  of 
Phobos  at  J2000.0,  N  is  the  motion  of  the  node  of  Phobos’ 
orbit,  and  iVis  the  position  of  the  node  of  Phobos  at  J2000.0. 
These  values  for  the  physical  quantities  for  Phobos  give  a 
precession  constant  of 


C 


Years  from  present 

Fig.  2.  Very  long-period  nutation  in  longitude  over  historical  periods 
of  time.  Over  short  time  spans  ±  10  000  yr  the  change  in  the  longitude 
of  Mars’  pole  is  nearly  quadratic  with  time.  The  coefficient  for  the 
quadratic  term  is  only  0"0001  yr  -  2 . 
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Q  =  —  =  ( 1,95  ±  0.15)  X  10~13  rad/s 

fir  q 

=  1"26±0"10  yr-1. 

Deimos: 

m  =  ( 1.8  + 0.15) X1015  kg, 
r  =  (2.3459  ±  0.0001 )  X  107  m, 
i  =  2!695  ±  0°030, 
n  =  1817.8545  ±  0.0003  rad/yr, 

A0  =  3.768  ±  0.001  rad, 

N=  -0.1156  +  0.0001  rad/yr, 

7V=  0.16  ±0.01  rad. 

These  values  give  the  precession  constant  due  to  Deimos  of 

3G^4=(2  i±02)x10-15  rad/s 
fir  q 

=  0"014  ±  07001  yr-1. 

Thus,  both  of  the  satellites  of  Mars  have  the  potential  of 
contributing  significantly  to  the  precession  and  nutation  of 
Mars. 

Since  both  of  the  satellites  orbit  Mars  in  planes  near  that  of 
the  planet’s  equator,  the  easiest  method  to  deal  with  them  is 
to  return  to  the  equations  of  motion  (4)  and  rewrite  them 

—  P  =  Q  sin  <5  cos  <5  cos  a  =  —  sin  26  cos  a, 


—  A  sin  e  —  Q  sin  S  cos  8  sin  a  =  —  sin  26  sin  a,  ( 16) 

2 


where  Q—  {hGmJ2) / (r'flq) .  Since  the  orbital  planes  of 
Mars’  moons  are  nearly  the  same  as  the  equatorial  plane  of 
Mars,  the  values  for  a  and  6  as  functions  of  time  are 


There  is  no  precession  term  contributed  by  the  Martian 
satellites  in  Eq.  (20)  because  the  invariable  planes  of  the 
satellites’  orbits  are  both  essentially  the  plane  of  Mars’  equa¬ 
tor.  The  similarity  of  the  two  planes  is  implicitly  included  in 
Eq.  ( 17)  for  6(t).  Thus,  the  precession  contributions  of  the 
satellites  are  cancelled  out  by  the  motions  of  the  nodes  of  the 
orbits. 

Evaluating  Eqs.  (19)  and  (20)  and  setting  the  constants 
of  integration  to  ^(0)  =0  and  e(0)  =  £0,  using  the  orbital 
parameters  for  Phobos  and  Deimos  gives 

Phobos: 

e  =  (0:'0040  ±  070003)  sin  [(  -  2°776  ±  0°003)/ 

+  (2°65  +  0707)  ]  +  (7"7±0"7X10~7) 

Xcos[  ( 14  395T12  ±  0.01  )t  +  (2!98  ±  0.07)]  +  e0, 
i/>  =  (474  +  0"2)  X 10“ 5  sin [(14  395!  12  +  0°01)/ 

+  (2798  ±  0707)]  +  (0"23  ±  0"01) 

Xcos[(  -  27776  +  0.003)/+  (2°65±0.07)]. 

Deimos: 

e=  (0"0028  ±0"0002)sin[(  -  0!1 156  ±  0°0001)/ 

+  (0°i6  ±  o°oi) ]  -  (9ro±o:'7xio-8) 

Xcos[  (3635°8264  ±  0.0004)/ 

+  (7.38  +  0.01 )  ]  +  £0, 

=  (1"9  ±  0'.'2X  10-6)sin[(3635!8264  ±  0.0004)/ 

+  (7738  +  0.01)] 

+  (0'.'060  ±  0.004)cos[ (  -0°1156±  0.0001)/ 

+  (0°16  ±  0.01)], 


a(t)xnt  +  / 10, 

6(t)  ~i  sin(«/  +  A0  —  Nt  —  N),  (17) 

where  i  is  the  inclination  of  the  satellite’s  orbit  to  Mars’ 
equator.  Substituting  these  relations  for  a  and  6  into  the  first 
of  Eqs.  ( 16)  and  assuming  /  is  small  gives  to  first  order: 

—  fix  ■Ysin2[/sin(«/  +  A0  —  Nt  —  JV)]cos(«/ +  A0) 

=  —  [sin(2n/  +  2A0  —  Nt  —  N)  +  cos  (Nt  +  N)  ]. 

2  (18) 

Integrating  Eq.  (18)  gives  the  nutation  in  latitude  as  a  func¬ 
tion  of  time  as 

6=  —  /?  =  —  (4-  sin(M  +  N) 

2  \N 

- - — t-  cos(2n/  +  2A0  —  Nt  —  N)  \  +  C.  (19) 

2n  —  N  J 

Performing  the  same  set  of  operations  on  the  second  of 
Eqs.  (16)  and  using  if>  =  —  A  gives  the  nutation  in  longitude 
as  a  function  of  time  as 

ifi  =  — (-  -1  v  sin(2n/  +  2A0  -  Nt  -  N) 

2  sin  6  \2n  —  N 

-4-sin(M  +  A0)  +  C.  (20) 

N  J 


where  the  time  /  is  in  years.  Only  one  term  in  each  of  these 
expressions  for  the  nutation  of  Mars  is  large  enough  to  be 
significant.  The  nutation  in  longitude  from  the  torque  of  the 
Martian  moons  is  deceptively  large.  The  reason  for  the  large 
nutation  in  longitude  is  because  it  is  proportional  to 
(sin  e)  “ 1 .  That  is,  the  motion  of  the  pole  in  longitude  moves 
a  large  distance,  but  along  a  very  small  circle  on  the  sky.  The 
main  source  of  uncertainty  for  the  nutation  from  the  satel¬ 
lites  is  the  inclination  of  Phobos’  orbital  plane  and  mass  for 
the  nutation  driven  by  Phobos  and  Deimos’  mass  for  the 
nutation  driven  by  Deimos.  Deimos,  despite  its  small  mass 
and  greater  distance  from  Mars  than  Phobos,  contributes  to 
the  nutation  because  of  the  slow  motion  of  its  node. 

The  third-order  term  for  the  expansion  of 
sin[2/sin(n/  +  /i0  —  Nt  —  Ar)]  in  Eqs.  (17)  is  at  least  3.7 
orders  of  magnitude  smaller  than  the  first-order  term  for  an 
inclination  of  >2°.  Thus,  the  higher-order  terms  in  the  ex¬ 
pansion  are  insignificant. 

The  nutation  that  arises  from  the  eccentricity  of  the  satel¬ 
lites’  orbits  is  all  insignificant  because  the  periods  for  the 
largest  nutation  components  are  on  the  order  of  1/n,  where  n 
is  the  mean  orbital  motion.  Since  the  amplitude  of  a  nutation 
is  inversely  proportional  to  its  period,  a  long-period  nutation 
has  a  larger  amplitude  than  a  short-period  one.  The  nutation 
components  arising  from  the  motion  of  the  nodes  are  signifi¬ 
cant  because  the  periods  are  long  (828  days  for  Phobos  and 
19  850  days  for  Deimos).  The  main  nutation  components 
arising  from  the  eccentricity  of  the  orbits  will  have  periods  of 
0.319  days  for  Phobos  and  1.26  days  for  Deimos.  Hence 
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these  nutation  components  will  be  about  the  same  size  as  the 
two  small  nutation  components  that  were  found  in  the  first 
order  expansion  in  Eqs.  (19)  and  (20),  and,  therefore, 
insignificant. 

5.  PRECESSION  CONTRIBUTIONS  FROM  THE  PLANETS 

Appendix  D  shows  that  the  perturbations  of  the  planets 
on  Mars’  orbit  cause  change  in  the  orbit  and  hence  a  change 
in  the  precession  rate  and  nutation  component  amplitudes 
from  the  change  in  the  gravitational  pull  of  the  Sun.  In  this 
section  the  direct  effect  of  the  gravitational  torques  of  the 
planets  on  Mars  will  be  assessed. 

Again  returning  to  Eq.  (4)  and  rewriting  its  first  two 
terms  gives  the  equation  of  motion  as 

3G  J i  ///  o  p 

—  p  — - -  —  sin  o  cos  o  cos  a 

ft  q  r 

—  Q  — l-  (2  sin  2 0  cos  e  +  sin  2A  sin  e 

+  sin  2/1  cos  2/1  .cos  e)  (21) 

and 

:  3  G  if,  ///  ■  c  o’ 

—  A  sin  e  = - -  — -  sin  o  cos  o  sin  a 

ft  q  P 

=  Q  — r  ( 3  cos  W  sin  2e  —  sin  2e 


—  cos  2/t  sin  2e  —  cos  2/3  cos  2A  sin  2e 


+  4  sin  2/3  sin  A  cos  2e),  (22) 


where  m  is  the  mass  of  the  perturbing  planet  and  r  is  the 
distance  between  the  planet  and  Mars.  Because  both  Mars 
and  the  perturbing  planet  are  circling  the  Sun  the  distance 
between  the  planet  and  Mars  varies  greatly  as  a  function  of 
time.  Thus,  the  distance  as  a  function  of  time  will  be  estimat¬ 
ed  using  a  series  approximation  in  Appendix  B. 

All  of  the  functions  on  the  right-hand  side  of  Eq.  (96)  in 
Appendix  E  have  the  common  factor 


2i  =  G 


m  sin  2c 
2ra>  sin  e 


(23) 


Substituting  in  the  appropriate  values  for  the  effect  of  the 
Earth  on  Mars  gives 

QlE  =  Ol'OOO  1044  — ^ - 1  000Ea  Sin50°40 

AU3  yr  2(1.524  AU)3  sin  25°20 

=  0"000 026  69  yr”1, 

where  Ea  is  the  mass  in  Earth  masses.  Substituting  in  the 
appropriate  constant  values  for  Jupiter  gives 

Qu  —  0"000  2132  yr-1. 


The  precession  in  longitude  is  determined  from  the  first 
full  term  on  the  right-hand  side  of  Eq.  (96): 


V  =  Q^B50/2(l+a2)t. 

4 


(24) 


Substituting  in  the  values  QIE,  a,  and  B  5()/2  for  the  torque  of 
the  Earth  on  Mars  gives 


=  (0°000  026  69  yr- l) (38.006) 
.000  AU) 


X 


[l  /LOOP  AUY]f  =  o„ooQ  363  yr 
L  \  1.524  AU/  J 


This  precession  is  almost  three  orders  of  magnitude  smaller 
than  the  uncertainty  in  the  solar  precession  alone  and,  hence, 
insignificant.  Substituting  in  the  similar  set  of  constants  for 
Jupiter  gives  a  precession  of 

Vj  =  0"000  198  yr  't. 

Again  the  precession  contributed  by  Jupiter  is  insignificant. 
Summing  the  contributions  from  all  of  the  planets  in  the 
solar  system,  the  total  precession  contributed  by  the  planets 
is  calculated  to  be 

yT  =  0"000  624  yr-  't. 

Thus,  the  Earth  and  Jupiter  contribute  89.9%  of  the  total 
planet  driven  component  of  Mars’  precession. 

Similarly,  the  nutation  in  longitude  components  all  have 
the  general  form  of 

Q^B5/2  i  ,  (25) 

4  n  sm  nt 

where  j  =  0,  1,  or  2,  and  n  is  the  sum  of  the  mean  motion  for 
each  of  the  perturbing  planets  and  Mars  times  some  integer 
coefficient.  Therefore,  based  on  the  size  of  the  coefficient 
Qie  the  amplitude  of  a  nutation  term  contributed  by  the 
Earth  will  be  significant  only  if 

B\/2~  >149.9. 
n 

With  the  correspondingly  larger  coefficient  Qu,  the  ampli¬ 
tude  of  a  nutation  term  contributed  by  Jupiter  will  be  signifi¬ 
cant  only  if 

B\/2—  >18.76. 
n 

The  Laplace  coefficients  for  both  planets  steadily  decrease 
with  increasing  k.  The  largest  coefficient  value  for  the  Earth 
is  B  q  2  =  38.006  and  the  largest  Laplace  coefficient  value  for 
Jupiter  is  B\/2  =  3.418.  The  value  of  a"  will  always  be  less 
than  or  equal  to  1.  Therefore,  if  a  nutation  term  contributed 
by  the  Earth  is  large  enough  to  be  significant  it  must  have  a 
mean  angular  motion  of 

n<  38.006/149.9  =  0.2535  rad/yr. 

And  a  significant  nutation  term  for  Jupiter  must  have 
n  <0.1822  rad/yr. 

For  the  Earth,  the  smallest  n  for  k  <  9  is  for  nE  —  3 
—  1.7511  rad/yr.  Therefore,  a  significant  nutation  where 
the  Earth  as  the  driving  body  does  not  exist.  However,  for 
Jupiter,  the  smallest  possible  n  for  k<  9  is 
3 ttj  —  n  #  =0.1 107.  Picking  out  the  related  terms  for  this 
value  of  n  gives  the  term 

q/y  =  Q^L^2±  J_  [(1  +a2)B5/2  +  aB5/2] 

2  r  sin  e  4 

X - - - sin(3  nj  —  n^)t.  (26) 

3  nj  -  n* 

Evaluating  this  equation  gives  a  nutation  term  of 
=  (0:'000  714  ±  0"000  002)sin(6°342t  +  315T1 1 1 ), 

where  t  is  in  years.  This  is  the  largest  nutation  in  longitude 
term  for  Mars  being  driven  by  Jupiter.  Although  it  is  a  factor 
of  3.2  larger  than  any  planetary  nutation  component  of  on 
the  Earth  (Vondrak  1982),  it  is  just  barely  significant. 

Like  Eq.  (96),  Eq.  (98)  in  Appendix  E  has  a  very  simple 
general  form  that  makes  the  search  for  nutation  in  latitude 
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components  that  are  worth  evaluating  simple.  The  general 
form  of  the  motion  in  latitude  terms  contains 

Q2Bo2  —  cosrcf,  (27) 

n 

where  Q2  —  Qm/2r\  sin  e.  Comparing  Q2  with  (?i  for  the 
nutation  in  longitude  equation  gives 

Qt/Q2  =  sin  2e/sin2  e.  (28) 

Substituting  in  the  present  obliquity  of  Mars  to  its  orbit  gives 

Qt/Q2  =  0.2658. 

Thus,  a  given  nutation  in  obliquity  term  is  only  about  27% 
the  size  of  the  equivalent  nutation  in  longitude  term.  There¬ 
fore,  the  one  nutation  in  obliquity  term  that  corresponds  to 
3 rij  —  n  y/  is  insignificant.  Its  amplitude  is 

8e=  (0"000  1897  ±  (TOOO  0005)  cos  (6)342t  +  315!111). 

The  direct  gravitational  attraction  of  the  planets  on  Mars 
does  not  contribute  significantly  to  the  precession  and  nuta¬ 
tion  of  Mars.  However,  as  shown  in  Sec.  3.3,  the  change  in 
the  orbit  of  Mars  arising  from  perturbations  by  the  planets 
does  cause  significant  changes  in  the  solar  precession  and 
nutation  of  Mars  over  very  long  periods  of  time. 

6.  CONCLUSIONS 

The  total  precession  in  longitude  for  Mars  is 

A>P=  —  7"296  ±  0"021  yr“'. 

This  precession  includes  two  components  resulting  from  ex¬ 
tremely  long-period  (>210  000  yr)  nutation  components 
from  long  period  changes  in  Mars’  orbit.  Aside  from  the  two 
contributions  from  the  change  in  Mars’  orbital  plane  with 
time,  the  precession  is  in  excellent  agreement  with  the 
precession  found  by  Reasenberg  &  King  ( 1979 ) .  These  same 
very  long-period  nutation  components  also  contribute  a 
“precession”  in  latitude  of 

Ae  =  0"4255  +  0"0012  yr”1. 

The  nutation  components  in  longitude  with  amplitudes 
larger  than  0"001  consist  of  eight  solar  terms  in  Table  2  and 
one  nutation  each  contributed  by  the  motions  of  the  nodes  of 
Phobos  and  Deimos.  The  nutation  components  in  obliquity 
consist  of  the  six  solar  terms  in  Table  2  and  one  nutation 
term  each  from  the  motions  of  the  nodes  of  Phobos  and  Dei¬ 
mos.  Again,  the  solar  nutation  in  both  longitude  and  latitude 
are  in  excellent  agreement  with  those  found  by  Reasenberg 
&  King  (1979). 

Although  the  torques  of  the  Earth  and  Jupiter  on  Mars  are 
several  times  larger  than  the  torques  of  the  planets  in  the 
solar  system  on  the  Earth,  only  a  single  nutation  in  longitude 
driven  by  Jupiter  is  significant  at  the  milliarcsecond  level. 

Looking  at  the  nutation  as  a  whole,  Fig.  3  shows  the  nuta¬ 
tion  of  the  Martian  pole  with  respect  to  the  mean  pole  over  a 
period  of  686.9  days  ( 1  Martian  sidereal  year).  This  figure  is 
very  different  from  the  same  one  for  the  Earth  (Fig.  4).  The 
difference  between  the  two  figures  is  the  result  of  the  differ¬ 
ence  in  the  main  driving  mechanisms  for  the  dominant  terms 
in  the  nutation  series  for  the  two  planets.  For  the  Earth,  each 
of  the  dominant  terms  is  a  factor  of  a  few  smaller  than  the 
next  larger  term  and  the  periods  are  different  by  a  factor  of 
several.  The  periods  and  amplitudes  are  the  result  of  the 
Earth  having  two  major  bodies,  the  Sun  and  the  Moon,  driv¬ 
ing  its  significant  nutation  components.  For  example,  the 
three  largest  nutation  components  in  longitude  for  the  Earth 


Nutation  in  Longitude  (") 


Fig.  3.  Path  of  Mars’  true  celestial  pole  with  respect  to  its  mean  celes¬ 
tial  pole.  The  path  of  Mars’  true  celestial  pole  on  the  sky  over  the 
period  of  1  Martian  sidereal  year  (686.93  days).  The  mean  pole  moves 
at  a  rate  of  —  7"3  along  the  small  circle  of  precession.  The  motion  of 
the  mean  pole  on  the  sky  is  about  3.1  arcsec  per  year.  The  “boomer¬ 
ang”  motion  of  the  true  pole  about  the  mean  pole  is  the  nutation. 


have  amplitudes  of  —  17"20,  —  1"32,  and  —  0"23  with  per¬ 
iods  of  6798,  183,  and  14  days,  respectively.  The  irrational 
ratios  of  the  periods  and  amplitudes  give  rise  to  the  ellipse  on 
ellipse  structure  for  the  Earth’s  nutation  in  Fig.  4.  For  Mars, 
however,  the  three  largest  nutation  components  in  longitude 
have  amplitudes  of  1"  10,  —  0"64,  and  —  0" 24  with  periods 
of  344,  687,  and  229  days,  respectively.  The  ratio  of  the  per¬ 
iods  is  nearly  2:1:3.  The  ratios  of  the  periods  are  nearly  in 
resonance  because  the  only  motions  involved  in  driving  the 
nutation  are  Mars’  orbital  motion  and  the  motion  of  the  peri¬ 
helion  with  respect  to  the  ascending  node  of  the  equator  of 
Mars.  The  motion  of  the  perihelion  with  respect  to  the 


Nutation  in  Longitude  (") 

Fig.  4.  Path  of  the  Earth’s  true  celestial  pole  with  respect  to  its  mean 
celestial  pole.  The  path  of  the  Earth’s  true  celestial  pole  in  the  sky,  over 
an  18  yr  period  beginning  at  J2000.0,  with  respect  to  the  mean  pole. 
The  mean  pole  moves  along  a  smooth  arc  at  a  rate  of  about  50"  per 
year  along  the  small  circle  described  by  precession.  The  motion  along 
the  small  circle  is  seen  as  a  motion  of  20  arcsec  per  year  on  the  sky.  The 
complex  epicyclic  motion  of  the  true  pole  is  nutation. 
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Table  3.  The  apparent  amplitude  of  the  combined  mean  motion  and 
perihelion  terms  for  the  three  largest  nutation  components  in  longitude. 


Period 

Amplitude 

Phase 

(at  J2000.0) 

(day) 

(") 

O 

686.82 

-0.5577 

263.27 

343.44 

1.131 

127.20 

228.97 

0.2433 

110.19 

equinox,  0°03  yr  " 1 ,  is  very  slow  in  comparison  to  the  orbital 
motion,  191°417  yr  1  (Bretagnon  1982).  The  very  slow 
motion  of  the  perihelion  causes  the  periods  of  the  dominant 
nutation  components  in  longitude  and  latitude  to  have  per¬ 
iods  that  are  nearly  simple  ratios.  This  simple  structure  is 
seen  on  the  sky  as  Fig.  3. 

The  slow  motion  of  the  perihelion  also  causes  the  nutation 
components  in  longitude  with  the  forms  of  sin  (2  A  +  nJf) 
to  appear  over  short  periods  of  time  (decades)  as  a  single 
nutation  with  periods  of  and  amplitudes  and  phases  of 

A  sin  (2A  +  nJf)  +  B  sin  nJ? 

=  JA2  +  B2  +  2 AB  sin  2A  cos 

x(nJf  +  tan~ 1  - C0S  +  — )  .  (29) 

At  J2000.0 
A  =  250°70. 

Table  3  gives  the  apparent  amplitudes  and  phases  for  the 
three  largest  nutation  components  in  longitude. 

APPENDIX  A:  THE  EQUATION  OF  MOTION 


The  first  four  Legendre  polynomials  used  in  the  equation 
of  motion  (2)  in  geometric  algebra  form  are 


/’o(rT')  =  1, 

(30) 

P,(  r-r')  =  r-r', 

(31) 

P2(  r-r')  =  ^[3(r-r')2  —  rV2], 

(32) 

P3  ( r-r' )  =  i  [  5  ( r-r' ) 3  -  3rV2  ( r V )  ] . 

(33) 

Integrating  Eq.  (2)  for  the  potential,  the  E^rr')  term 
disappears.  The  gravitational  potential  for  a  planet  is  then 

(34) 

where 

S(r)  =  r2  Tr(/)  —  3r-/(r), 

(35) 

where  1  is  the  inertia  tensor  and  Tr(/)  is  the  trace  of  the 
inertia  tensor. 

The  gravitational  field  of  an  extended  body  is  found  by 
taking  the  gradient  of  Eq.  (34) 

g(r)  =  -  V«f>(r) 

— +■••)]■ 

=  (uW  -  ^  [5(f)  -  ff-5(f)f  ]  +  •••). 

(36) 

Substituting  for  S(r)  the  gravitational  field  becomes 


g(r)  = 


-GJ( 

r2 


r  —  3J2(yJ  { [ 5 ( r-e ) 2  —  l]r 


2[f-e]e}  + 


(37) 


where  R  is  the  equatorial  radius  for  the  body,  e  is  the  vector 
for  the  axis  of  the  principal  moment  of  inertia  about  the  pole 
of  the  assumed  hydrostatic  spheroid  for  the  planet,  r  is  the 
unit  vector  from  the  planet  to  the  torquing  body,  and  J2  is  the 
dynamical  form  factor  defined  by 


(38) 


where  /3  is  the  principal  moment  of  inertia  about  the  axis  of 
figure  and  /  is  the  equatorial  principal  moment  of  inertia  for 
a  spheroidal  planet. 

Thus,  the  force  of  a  point  mass  m  on  an  extended  mass  , 
is 

r-  3/2^-j  {[5(r-e)2-  l]r 


F  =  ^g(r)  = 


—  m 

? 


-  2[r-e]e}  +  •■•  j  . 


(39) 


The  torque  of  m  on  J/  is  perpendicular  to  the  vector  from 
^  to  m,  that  is  r  =  r  X  F,  so  all  of  the  terms  in  the  direction 
of  r  can  be  ignored.  The  first-order  perpendicular  compo¬ 
nent  of  F  is  simply 

Fis  — 3^^j!(K)e.  (40) 


The  vector  for  the  torque  per  unit  moment  of  inertia  H  is 

H  =  —  =s  -3-^L^(f*e)e.  (41) 

I  r3  q 


The  steady  precession  for  a  torquing  body  whose  position 
is  fixed  with  respect  to  the  planet  is 


ft, 


- H 
fte 


(42) 


Substituting  in  Eq.  (41 )  for  the  torque  per  unit  moment  of 
inertia  gives 


ft, 


3  Gm  J-,  r-e  „ 

- ; — - e. 

ftr3  q  r 


(43) 


This  motion  is  the  mean  precession  of  Mars  assuming  that 
the  position  of  the  body  providing  the  torque  is  fixed  with 
respect  to  the  planet.  Mars,  however,  is  moving  with  respect 
to  the  other  masses  in  the  solar  system,  so  the  time  averaged 
value  of  (r-e)/r  is  needed.  The  vector  r  is  a  function  of  time 
so  it  is  now  convenient  to  introduce  coordinates  that  will 
represent  r(t)  and  can  be  used  to  produce  measurable  quan¬ 
tities  for  the  precession  and  nutation  of  Mars.  The  coordi¬ 
nate  system  most  useful  for  determining  [r(t)-e]  and  the 
coordinate  system  most  suited  for  measuring  precession  and 
nutation  are  rotated  with  respect  to  one  another.  First,  the 
celestial  coordinate  system  for  the  position  of  a  perturbing 
body  with  respect  to  Mars  is  the  equatorial  coordinate  sys¬ 
tem,  the  Martian  Right  Ascension  a,  and  Martian  declina¬ 
tion  S  of  the  perturbing  body.  These  coordinates  are  defined 
to  be  exactly  analogous  to  the  Earth’s  equatorial  celestial 
coordinate  system.  The  epoch  used  for  the  coordinate  sys¬ 
tem  is  J2000.0.  In  this  system  of  coordinates  the  dot  product 
(r-e)/r  is  simply 
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r(t)-e/r  =  sin  <5(0-  (44) 

The  direction  vector  e  becomes 
e  =  cos  8(t)a,  (45) 

where  a  is  the  direction  vector  in  curvilinear  coordinates  of 
increasing  a.  Finally,  since  the  torque  is  the  cross  product  of 
F  and  r,  the  motion  is  in  the  a  direction.  Substituting  these 
values  for  (r-e)/r  and  e  into  Eq.  (43)  gives 

H  =  — ’Ll  sin  S(t)  cos  S(t)a.  (46) 

r  q 

The  second  coordinate  system  used  is  for  the  measure¬ 
ment  of  precession  and  nutation.  Precession  and  nutation 
involve  a  change  in  the  orientation  of  the  rotation  axis  with 
time,  so  it  is  preferable  to  describe  their  components  in  a 
coordinate  system  that  is  independent  of  the  position  of  the 

I 


planetary  pole.  The  coordinate  system  usually  used  is  the 
Euler  angles,  tp,  <f>,  and  e.  The  angle  ip  is  the  angle  from  the 
equinox  of  a  chosen  epoch  to  the  intersection  of  the  ascend¬ 
ing  node  of  Mars’  equator  of  date  along  the  plane  of  the  orbit. 
The  angle  <j>  is  the  rotation  of  a  point  on  the  surface  of  Mars 
about  the  axis  of  rotation.  The  angle  e  is  the  inclination  of  the 
equator  of  Mars  to  the  plane  of  its  orbit  about  the  Sun. 

Using  the  Euler  angles  and  the  ecliptic  coordinates  the 
angular  momentum  vector  of  the  planet  is 

1  =  /3fl[  (cos  P  cos  A), (cos  P  sin  A  cos  e  —  sin  P  sin  e), 

(cos  P  sin  A  sin  e  +  sin  P  cos  e)],  (47) 

where  A  is  the  ecliptic  longitude  of  the  pole  of  rotation  and  P 
is  the  ecliptic  latitude  of  the  pole.  Thus,  the  time  derivative  of 
the  angular  momentum  in  the  ecliptic  coordinate  system  is 


1  =  /3fl{  [  —  p  sin  p  cos  A  —  A  cos  p  sin  A  ] , 

[/?(sin/3sin  A  cos  e  +  cos  P  sin  e)  +  A  (cos  P  cos  A  cos  e)  —  e(cos  p  sin  A  sin  e  +  sin  P  cos  e  )  ], 

[/(cos  P  cos  e  —  sin  P  sin  A  sin  e)  -f  A ( cos  P  cos  A  sin  e)  -f  e(cos/J  sin  A  cos  e  —  sin  p  sin  e)  ]}  .  (48) 


If  the  epoch  of  the  ecliptic  coordinate  system  is  chosen  to  be 
near  the  time  of  observation,  then  the  angular  momentum 
vector  will  point  near  A  =  90°,  p  =  90°  —  e,  and  the  obliqui¬ 
ty  e  is  approximately  constant  so 

cosAsO,  sinA~l, 
cos  P  ~  sin  e,  sin  p  ~  cos  e, 
and 

exO.  (49) 

Using  approximations  (49)  in  Eq.  (48)  the  time  derivative 
of  the  angular  momentum  vector  becomes 

iss/3fl[  —  A  sin  e,  —  j3( cos2  e  +  sin2  e), 

/(sin  e  cos  e  —  cos  e  sine)  ] 

=  /3n(  -A  sine,  -/,0).  (50) 

It  is  now  possible  to  resolve  H,  Eq.  (43),  into  its  compo¬ 
nent  parts  and  equate  it  with  1.  Since  1  =  F  =  JAW 

(  —  A  sin  e,  —  /, 0) 

_  3 Gm  ^2(sjn  g  CQS  g  sjn  a  sjn  g  cos  g  CQS  ao).  (5 1 ) 

r  Cl  q 

The  equatorial  coordinates  used  for  the  position  of  the  body 
supplying  the  torque  needs  to  be  transformed  to  the  ecliptic 
coordinate  system.  The  transformation  between  the  two  co¬ 
ordinate  systems  is  given  in  The  Explanatory  Supplement  to 
The  American  Ephemeris  and  Nautical  Almanac  ( 1974) . 


APPENDIX  B:  EXPANSION  OF  THE  EQUATION  OF  MOTION 
FOR  AN  ELLIPTICAL  ORBIT 

The  process  of  expansion  for  an  elliptical  orbit  is  begun  by 
dropping  the  assumption  of  uniform  motion  and  distance 
from  the  attracting  body.  Equations  (6)  and  (7)  are  rewrit¬ 
ten 


r 


A  sin  e  —  ®  sin  e  cos  e  sin 2(f  —  A), 

r 

P  =  Q  sin  esin(/—  A)cos(/—  A),  (52) 

r 


where  Q  =  ( 3  /77.GJ2) / (Wq)  and  f  is  the  true  anomaly  of 
Mars.  Using  standard  trigonometric  identities  the  sine  and 
cosine  functions  in  Eq.  (52)  are  rewritten 

sin2(/  —  A)  =  \  —  \  cos  2/cos  2A  —  |  sin  2/sin  2A 
and  (53) 

sin(/—  A)cos(/—  A)  =  1  sin  2/cos  2A  —  \  cos  2/sin  2A. 

Neither  the  true  anomaly  nor  the  distance  from  the  Sun,  r, 
are  linear  functions  of  time.  However,  both  the  true  anomaly 
and  the  distance  r  can  be  represented  by  series  expansions, 
Smart  (1953)  Chap.  3: 

3  =  X  Dn  cos (nJT), 

'  n  =  0 

^jsin2/=  jr  En  sin(nJ/),  (54) 

r  n  =  0 

^-cos2/=  ^  Fn  cos(nJf), 

r  n  =  0 


where  a  is  the  mean  distance  of  the  planet  from  the  Sun.  The 
coefficients  for  these  series  have  been  tabulated  by  Jarnagin 
( 1965)  and  the  series  for  Mars,  truncated  to  six  significant 
figures  in  the  eccentricity  e,  are  given  in  Table  4.  The  value 
used  for  the  mean  eccentricity  e  is 

<?  =  0.093  400  619  9474  (Bretagnon  1982). 


Substituting  Eqs.  (54)  into  Eqs.  (53)  and  then  substitut¬ 
ing  that  result  into  equations  of  motion  ( 52 )  and  simplifying 
gives  the  motion  in  longitude  and  latitude  as 


3  - Q 

A  =  — cos  e 


—  +  X  (~ —  cos  2A^ 

2  „  =  ,l  2  2  ) 


B  \ 

Xcos (nJf) - sin  2A  sin (nJf)J  , 
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Table  4.  Coefficient  series  in  eccentricity,  e ,  for  polynomial  expansions 
in  the  mean  longitude,  . 


Series 

Trig.  Func. 
of  L 

Coefficient  Series 

Value 

(a/r)3 

cos  2f 

cos  1 L 

-e  +  3/48  e3  -  5/384  e5 

-0.0466495 

cos  2L 

1  -  5/2  e2  +  13/16  e4  - 

35/288  e6 

0.9782526 

cos  3  L 

7/2  e- 123/16  e3  + 

489/128  e5 

0.3206656 

cos  4  L 

17/2  e2- 115/6  c4  + 

601/48  e6-  1423/360  e8 

0.0727009 

cos  5  L 

845/48  e3  -  32525/768  c5  + 

1457575/43008  e7 

0.0140167 

cos  6  L 

533/16  e4-  13827/160  e6  + 

728889/8760  e8 

0.0024783 

cos  1L 

228347/3840  c5  - 

3071075/18432  e7 

0.0004123 

cos  8  L 

73369/720  e6- 

775727/2520  e8 

0.0000659 

cos  9L 

12144273/71680  e7 

0.0000105 

cos  10  L 

634085  e8 

0.0000016 

(a/r)3 

sin  2f 

sin  1 L 

1/48  e3+  11/768  e5 

0.0000171 

sin  2  L 

2/48  e4  +  42/1440  e6 

0.0000032 

sin  3  L 

81/1280  e5 

0.0000004 

sin  4  L 

647720  c6 

0.0000001 

(a/r)3 

cos  0  L 

1  +  3/2  e2  +  15/8  e4  + 

35/16  e6 

1.0132297 

cos  1 L 

3  e  +  27/8  e3  +  261/64  e5 

0.2829808 

cos  2  L 

9/2  e2  +  7/2  e4  + 

141/32  e6 

0.0395258 

cos  3  L 

53/8  e3  +  393/128  e5  + 

24753/51 20  e7 

0.0054202 

cos  4  L 

77/8  e4  +  l  29/80  e6 

0.0007336 

cos  5  L 

1773/128  e5 

0.0000985 

cos  6  L 

3167/160  e6 

0.0000131 

cos  1L 

432091/15360  e7 

0.0000017 

E  =  -^sine/  jr  -^-cos2Asin(«.Z’) 

Cl  \n  =  0  2 

—  ^y-sin  2A  cos(n^)j  .  (55) 

The  position  of  the  perihelion  with  respect  to  the  Martian 
vernal  equinox  is  the  difference  between  the  advance  of  the 
perihelion  and  the  precession  of  Mars.  Using  the  Bretagnon 
( 1982 )  mean  orbit  for  Mars  and  the  circular  orbit  precession 
previously  determined,  the  position  of  the  perihelion  is 

A  =  94"050  -  (  -  7"488)  =  1017538  yr"1. 

The  motion  of  precession  also  folds  in  the  effect  of  the  mo¬ 
tion  of  Mars’  perihelion  with  respect  to  its  vernal  equinox. 
However,  the  change  in  the  precession  rate  from  the  circular 
orbit  value  is  too  small  to  affect  the  rate  of  precession  from  its 
inclusion  in  A  or  the  amplitudes  or  periods  of  the  nutation 
components. 

Collecting  like  terms  together,  constant  Eqs.  (55)  inte¬ 
grate  to  give  the  nutation  in  longitude  and  latitude  as 

A  =  cos  e(A0t  4-  ]T  -L-  [A„  sin  (nJf)  ] 

Cl  \  n  —  \  H 

B  \ 

-  C„  cos  (2A  +  nJf)  +  — sin(2A  +  n&)  , 


/?  = — sin  el  f  ——7—  cos  (2A  +  nJf) 
a  \,,  =  o  nJf 

+  -~sin  (2A  +  .  (56) 

The  coefficients  A0,  An,  B„,  and  C„  are  linear  combinations 
of  the  coefficients  in  Eq.  (54). 


APPENDIX  C:  PRECESSION  AND  NUTATION  FROM  CHANGES 
IN  ECCENTRICITY 

The  Bretagnon  ( 1982)  theory  for  Mars  uses  the  canonical 
elements  k  =  e  cos  Ft  and  h  =  e  sin  II,  where  e  is  the  eccen¬ 
tricity  and  n  is  the  longitude  of  perihelion.  Both  of  the  ele¬ 
ments  k  and  h  are  given  in  the  form 

k  =  '£TA"(Bn  sin  y„  +  Cn  cos  yn ), 


h  =  '2JDn{En  sinr„  +F„  cosy„),  (57) 

rt 

where  T  is  the  time  in  thousands  of  Julian  years  from 
J2000.0,  An  and  Dn  are  the  integer  exponential  coefficients 
of  the  time  dependence,  Bn,Cn,En,  and  Fn  are  the  amplitude 
coefficients  for  the  sine  and  cosine  functions,  and 

Yn  =  KnXZ  +  Kn2  V  +  KniE  +  KnttM  +  K„5J 

+  Kn6S  +  Knl  U  +  Kn%N,  (58) 

where  the  variables  Z,  V,  E,  M,  J,  S,  U,  and  A'' are  the  mean 
ecliptic  longitudes  of  Mercury,  Venus,  Earth,  Mars,  Jupiter, 
Saturn,  Uranus,  and  Neptune,  respectively,  and  the  Knm  are 
integer  coefficients  of  the  mean  longitudes.  Only  the  mean 
ecliptic  longitudes  of  the  Earth,  Mars,  and  Jupiter  are  in¬ 
volved  in  those  terms  which  may  be  large  enough  to  signifi¬ 
cantly  affect  the  amplitudes  of  Mars’  nutation  components. 

For  a  given  term  n  in  the  series  making  up  the  periodic 
terms  in  eccentricity,  the  conversion  from  h„  and  kn  to  e\  is 

e'  =h2n  +k2„  =G„  +H„  cos(2r„  +  *„),  (59) 

where 

Gn=\(Bl+E2n-Cl-Fl), 

Hn  =  ](C2„  +  F\  -  B\  -  El  )2  +  (B„Cn  +  EnFn  )2, 

and 


Xn 


=  arctan 


2(J?„C„  +EnFn) 

Bl+El-Cl-Fl  ‘ 


(60) 


The  conversion  from  h  and  k  to  e2  has  a  constant  offset  in 
eccentricity  and  a  cosine  function  with  an  angular  motion 
that  is  twice  that  of  the  original  angular  motion  in  h  and  k. 
The  eccentricity  is  determined  from  e2  by  applying  a  Taylor 
series  expansion  to  Eq.  (60).  To  first  order  the  expansion  is 


e=(Gn  +Hn  cos  j„) 1/2 


Hn  sin  Xn _ 

2(G„  +Hn  cos Xn)'/2Vn 

(61) 


In  Bretagnon’s  theory  for  Mars  there  are  seven  terms  in  h 
and  k  with  A„  =  D„  =  0  for  which  the  zeroth-order  Taylor 
series  term  for  eccentricity  is  larger  than  0.000  01  and  there 
is  one  term,  with  An  —  Dn  =  1,  for  which  the  zeroth-order 
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term  is  larger  than  0.000  001.  Inserting  the  values  of  these 
terms  into  Eq.  (41 )  shows  that  the  largest  first-order  term  in 
Eq.  (41 )  is  a  factor  of  1 1 .4  smaller  than  the  smallest  zeroth- 
order  term.  Hence  only  the  zeroth-order  terms  are 
significant. 

APPENDIX  D:  PRECESSION  AND  NUTATION  FROM  CHANGES 
IN  THE  ORBITAL  PLANE 


Because  the  change  in  obliquity  is  caused  by  changes  in 
Mars’  orbital  plane,  a  constant  plane  needs  to  be  chosen  to 
describe  its  motion.  The  inertial  reference  plane  chosen  is  the 
ecliptic  plane  of  the  Earth  at  J2000.0  and  the  variables  deter¬ 
mining  the  orientation  of  Mars’  orbital  plane  are  the  inclina¬ 
tion  of  the  orbit  to  the  ecliptic  i  and  the  position  of  the  as¬ 
cending  node  of  the  orbit  on  the  ecliptic  ft.  The  general 
development  of  the  long-period  motion  of  the  pole  will  fol¬ 
low  that  of  Ward  (1974)  and  Reasenberg  &  King  (1979). 
The  main  results  of  this  section  are  ( 1 )  a  “linear”  term  add¬ 
ed  to  the  precession  as  a  result  of  the  changes  in  the  orbital 
plane  and  ( 2 )  a  “linear”  term  for  the  change  of  the  obliquity 
in  time.  Both  of  these  terms  are  the  first-order  approxima¬ 
tions  of  sinusoidal  motions  which  are  approximately  linear 
over  periods  on  the  order  of  10  000  yr. 

The  first-order  equations  of  motion  are 


=  sin  i  cos  i/y/a  cos  e0  +  )  +  sin  ip0  — , 

dt  dt  )  Y0dt 

d\b,  dCl  .  . 

— —  =  —  at  cos  e, - sin  i  cot  en 

dt  dt  0 

X  sin  ipja  cos  e0  +  4-  cot  e0  cos  ipu  — 

\  dt  /  dt 


(62) 


where  e0  and  ip0  are  the  zeroth-order  solutions  for  the  motion 
e  =  e0 

and  (63) 

if>0  —  at  cos  e0  +  ft0  —  ft(f), 

where  'P  and  ft0  are  the  values  of  ^  and  ft  at  some  given 
initial  time  t0  and  ft  ( 0  is  the  value  of  ft  as  a  function  of  time. 
The  motion  in  latitude  is  rewritten 


^4  =  a  cos  e0[p  cos  (at  cos  e0  —  ft0  —  'P) 
dt 

—  q  sin(af  cos  e()  —  ft0  —  'P )  ]  +  —  sin  i  sin  Tp0, 

dt 


(64) 


where p  =  sin  i  cos  ft(f)  and  q  —  sin  i  sin  ft  (7),  the  canoni¬ 
cal  elements  (Bretagnon  1982)  for  the  orientation  of  the 
Mars’  orbit.  This  allows  the  motion  in  latitude  to  be  solved 
for  as  a  function  of  the  change  in  the  orientation  of  Mars’ 
orbit  with  time. 

The  canonical  elements  q  and  p  are  also  functions  of  time. 
The  first-order  solutions  for  p  and  q  are  represented  by 

<lf='LAJf  sin(®/  +  <5,  )> 

j 

Pf  =  cos ( +  SJ ) .  (65) 

j 


The  constants  Ag,  ttfj,  and  Sj  have  been  solved  for  the  major 
solar  system  bodies,  excluding  Pluto,  by  Brouwer  &  van 
Woerkom  (1950).  These  constants  for  Mars,  brought  for¬ 
ward  from  the  ecliptic  and  equinox  of  B1950.0  to  the  ecliptic 
and  equinox  of  J2000.0  and  modified  for  the  IAU  1984  val¬ 
ues  for  the  reciprocal  masses  of  the  planets  (Kaplan  1981), 
are  given  in  Table  5. 

Substituting  Eqs.  (65)  for  p  and  q  into  Eqs.  (64)  gives 
=  a  cos  e0^A4J  cos  ( b}  t  +  at  cos  e0 

dt  j 

+  Sj  —  ft0  —  *P)  +  —  sin  i  sin  if>0.  (66) 

dt 

Equation  (66)  is  integrated  to  obtain  the  first-order  equa¬ 
tion  for  the  nutation  in  latitude  of  Mars  as  a  function  of  time 
resulting  from  the  change  in  the  orientation  of  Mars’  orbit 
with  time. 


e,  =  a  cos  e0V - — - 

j  a  +  a  cos  e0 

X  sin (atjt  +  at  cos  e(l  +  Sj  —  ft0  —  T) 

+  sin  / sin  if>0  +  C,  (67) 

where  C  is  the  constant  of  integration.  Since  the  obliquity  to 
the  ecliptic  at  t  =  0  is  e0,  the  value  of  the  constant  of  integra¬ 
tion  is 

C=e0-~y - — ~ ^4j - sin  (<5  —  'P  —  ft0)  (68) 

j  fit  j  +  a  cos  e0 


and 


_  v-  Oj-*v 

ei  —  eo  —  X  ; 

j  Bj  +  a  cos  e0 

X  [sin(5T;t  -f  at  cos  e0  +  Sj :  —  'P  —  ft0) 

+  sin(<5;  —  T  —  ft0) ]  .  (69) 

Plotting  this  equation  for  ?!  in  Fig.  5  shows  a  change  in 
obliquity  with  time  that  appears  to  be  a  beat  frequency  of  two 
sine  waves  with  periods  of  1.16X  105  yr  and  1.29X  105  yr. 
These  are  the  frequencies  of  the  perturbations  of  Jupiter  and 
the  Earth  which  supply  83.1%  of  the  amplitude  to  the  per¬ 
turbation  in  latitude.  Over  short  periods  of  time  near  the 
present  era,  +  10  000  yr  (Fig.  1 ),  the  motion  in  obliquity  of 
the  Martian  pole  as  a  result  of  the  change  in  orientation  of 
Mars’  orbit  can  be  represented  by  the  linear  function 


Table  5.  Constants  for  the  canonical  elements  for  Mars. 


Planet 

(i) 

ni 

(”/vr.) 

8i 

<°) 

A4i 

(rad.) 

A4i 

(") 

i 

-5.202 

271.99 

0.0017940 

370.04 

2 

-6.571 

209.97 

0.0017989 

371.05 

3 

-18.744 

147.13 

-0.0359444 

-7414.06 

5 

-17.633 

188.68 

0.0502514 

10365.1 

6 

-25.734 

19.22 

0.0096568 

1991.86 

7 

-2.903 

207.44 

-0.0012561 

-259.09 

8 

-0.678 

95.00 

-0.0012286 

-253.42 

Planets  are  numbered  by  their  increasing  mean  distance  from  the  Sun 


©  American  Astronomical  Society  •  Provided  by  the  NASA  Astrophysics  Data  System 


■o 

liH 

,LO 


1522 


■CM 

IO 


i1^) 

ICTi 
I  CT» 


1522  JAMES  L.  HILTON:  MARS’  POLE 


Years  from  present 


7"  1 57  X 10“  6  yr-2  is  insignificant.  The  uncertainty  in  the 
linear  term  which  appears  as  a  “precession”  in  inclination  in 
the  final  solar  nutation  series  is  +  0'.'0012.  The  main  source 
of  uncertainty  is  the  value  of  the  inertia  ratio. 

Substituting  Eq.  (69)  for  e,  into  Eqs.  (62)  the  first-order 
equation  of  motion  in  longitude  becomes 


#i 

dt 


Bj  +  a  cos  e0 


X  [sin(5r/  +  at  cos  e0  +  S}  -  -  fl0) 


Fig.  5.  Very  long-period  nutation  in  obliquity.  The  change  in  the  incli¬ 
nation  of  Mars’  pole  with  respect  to  its  orbital  plane  resulting  from  the 
change  in  Mars’  orbital  plane  with  time. 


e  =  (25°2  ±  0!01 )  +  0"4255t,  (70) 

where  t  is  the  time,  in  years,  from  J2000.0.  This  linear  equa¬ 
tion  is  obtained  by  doing  a  quadratic  fit  to  function  ( 69 )  for  e 
over  the  period  J2000.0  +  10  000  yr.  The  quadratic  term  of 


+  sin  (<5,  —  'I'  —  fl0)  ]  1  — - - sin  /  cot  e0 


X  sin  ifiJa  cos  e  +  -^Q. )  +  cot  e0  cos  ifi0  — .  (71 ) 
v  dt  J  dt 

Approximating  cos  i~  1,  putting  all  of  the  time  derivatives 
on  the  left  side  of  the  equation,  and  substituting  Eq.  ( 63 )  for 
if>0  into  Eq.  (71)  gives 


—  +  Cl  —  sin  i  cot  e0  cos  ip0) 

dt 

(  _  B:N:  \ 

=  —a  cos  e0  —  V - — - [  sm(Bt  +  at  cos  e0  +  <5  -  -  fl0)  +  sin  ( <5,  —  -  fl0)l  I 

V  j  Bj  +  a  cos  e0  V 

+  a  sin  / cot  e0  cos  e0  sin  [fl(l)  +  at  cos  e0  —  T  —  ft0]  .  (72) 

The  above  equation,  like  Eq.  (64)  for  e„  depends  on  (sin  i  sin  fl),  so  the  same  set  of  identities  for  p  and  q,  Eqs.  (65),  are 
substituted  to  give 

+  fl  —  sin  /  cot  e0  cos  ^0)  =  —  a  cos  e%  cos^  — ^  sin(/q  -f  /z,  )^  —  a  sin  sin^  sin(/cy  +  fij )  j 

+  a  sin  /  cot  e0  cos  e0  ^Nj  sin (Bjt  +  at  cos  e„  +  <5,  —  T  —  fl0)  ,  (73) 


where 

e*  =  e0  +  X  ~J~J~  Sln/2 3f  > 

J  kj 

Kj  —  Bj  +  a  cos  e0, 
and 

fXj=Sj-^-  ci0. 

The  sine  and  cosine  functions  in  Eq .( 7  3 )  are  approximat¬ 
ed  in  Ward  ( 1974)  by 

b,N ,  .  \ 

C0S(X  Sln(  V  +  )  J  =  1  (75) 

and 

sin(  X  sin(«v^  +  =  X  ~L~L  sin(  V  +  )  • 

(76) 

However,  the  total  amplitude  for  the  second-order  terms  for 
the  cosine  function  approximation  is 


BjNj 


=  0.015  958  143 


2  i 

which  when  multiplied  by  a  cos  e*  gives  a  total  amplitude  of 
(T2422  +  0"0007. 


Similarly,  the  total  amplitude  for  the  third-order  term  in  the 
approximation  for  the  sine  function  is 


BjNj 


3 

=  0.000  950  3144. 


Multiplying  this  by  a  sin  e*,  the  total  amplitude  of  the  third- 
order  term  is 


0"003  393  ±  0"000  009. 

Thus,  both  the  second-  and  third-order  terms  of  the  sine  and 
cosine  functions  are  potentially  large  enough  to  make  a  sig¬ 
nificant  contribution  to  the  precession  of  Mars.  Therefore, 
the  equation  of  motion  in  longitude  resulting  from  perturba¬ 
tions  in  Mars’  orbit  will  be  approximated  to  third  order  by 
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—  ( if/,  +  fi  —  sin  /'  cot  e0  cos  if>0)  =  —  a  cos  £$ 
dt 
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-  a  sin  eff  Wj^j  sin (Kjt  +  fij)  +  X  ^  Bj^J  sin (Kjt  +  J 

+  a  cos  e0  cot  e0^Nj  sin {Kjt  +  fij ) . 


(77) 


D.  1  Zeroth  and  First-Order  Longitude  Terms 

Equation  ( 77 )  is  the  same  as  that  derived  by  Ward  (1974) 
with  the  addition  of  the  second-  and  third-order  terms  in  the 
approximation  for  the  sine  and  cosine  terms.  The  second- 
and  third-order  terms  in  Eq.  (77)  are  related  to  the  zeroth 
and  first-order  terms  by  addition,  so  they  will  be  integrated 
separately  and  the  results  added  to  the  first-order  terms. 

Keeping  the  zeroth-  and  first-order  terms  only,  Eq.  (77) 
becomes 

—  {ipx i  +  ft  —  sin  /  cot  e0  cos  if>0) 
dt 

—  —  a  cos  e*  —  a  sin  e*^  — ^ -  sin  («v,  r  +  fij ) 

j  Kj 


+  a  cos  e0  cot  O^Nj  sin  {Kjt  +  fij ) . 


(78) 


Integrating  Eq.  (78)  gives 
ifin  f  fi(f)  —  sin  i  cot  e0  cos  ifi0 


®,N, 


=  —  at  cos  ej  +  a  sin  e*Y  — - — —  cos  {Kjt  +  fi, ) 

J 


N, 


—  a  cos  en  cot  e0^  — —  cos  {Kjt  +  fij)  +  C.  (79) 


i  nj 

Evaluating  Eq.  (64)  for  e*  gives 
e*  =  25!  199  +  0°02. 


This  value  is  indistinguishable  from  e0.  The  values  of  the 
other  angles  at  Z(0)  are  fl(0)  =  fi0  and  ^,,(0)  =  Thus, 
the  constant  of  integration  is 

O -N  ■ 

C'  =  —  V — 1  1  («  +  sjj  cos  e0 )  cos  fij.  (80) 

.  j  sin 

The  first-order  approximation  for  t/rl{t)  is 


I - 

tpu  —  'M  —  at  cos  ej  +  fi0  —  fi (r) 

+  y  .  ,  [  (a  +  Bfj  COS  e0) 

j  sin  erf 

XCOS(Kjt  +  fij)  +  COS  Hj]  .  (81) 

The  second  term  on  the  right-hand  side  of  Eq.  (81)  is  the 
already  determined  precession  of  Mars.  The  fourth  term  is 
the  motion  of  the  node  of  Mars’  orbit  along  the  ecliptic  of 
J2000.0.  The  first  and  third  terms  in  Eq.  (81 )  are  the  posi¬ 
tions  of  the  ascending  node  of  Mars’  orbit  on  its  equator  and 
the  node  of  the  orbit  of  Mars  at  J2000.0  with  respect  to  the 
ecliptic  and  equinox  of  the  Earth  at  J2000.0. 

The  last  term  on  the  right  side  of  the  equation  is  the  first- 
order  motion  of  the  longitude  caused  by  the  changes  in  the 
orientation  of  Mars’  orbit  with  time.  As  with  the  motion  in 
latitude,  the  motion  in  longitude  is  represented  as  a  series  of 
nutation  components  with  periods  varying  from  7 1  42 1  yr  to 
1270  000  yr  (Fig.  6).  The  amplitudes  of  these  changes  in 
longitude  are  evaluated  to  range  from  32  200"  from  the  per¬ 
turbations  by  Jupiter  to  65"7  from  the  perturbations  by  Nep¬ 
tune;  however,  88.5%  of  the  total  amplitude  of  these  pertur¬ 
bations  are  the  result  of  the  influences  of  Jupiter,  the  Earth, 
and  Venus.  Over  short  periods  of  time  (  <  10  000  yr)  the 
zeroth-  and  first-order  parts  of  the  first-order  solution  for 
precession  and  nutation  in  longitude  appear  as  a  contribu¬ 
tion  to  the  precession  of 

i[>n  =  (0127  ±  0”0006  yr"1 
shown  in  Fig.  2. 

D.  2  Second  and  Third-Order  Longitude  Terms 

Carrying  out  the  squaring  for  the  second-order  cosine 
term  in  Eq.  (77)  and  simplifying  it  gives 


a  cos  e* 

—  a  cos  e* 


+ 


1(2 

y  jx  y  (yy~)  [ 1  - cos  2(V  +  2Vj  >  ] 

XX  ^cos  [ +  K>‘)t + y-i  +M  - c°s[(*)  -Kk)t+ Hi  - Mk ] , 


(82) 


where  k  runs  from  j  +  1  to  8.  This  is  then  integrated  to  obtain 
a  cos  e* 

2 


r  [?{t  [(^)  -^-sin  2('',+ HI 

+  XXr^A  [ 1 sin  [(*■;  +Kk)t  +  fij  +fik] - l- sia[(Kj  -  Kk)t  +  fij  -  fik]  ) 

j  k  \  KjKk  (  Kj  -j-  Kk  Kj  —  Kk  J  / 


+  c. 


(83) 
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The  first  term  in  Eq.  ( 83 )  evaluates  to  be  a  secular  precession  of 
8xl>X2  =  0"0797  +  or0OO2  yr”1. 

The  other  three  terms  in  Eq.  (83)  are  63  periodic  components  with  amplitudes  as  large  as  0"0650  with  periods  ranging  from 
44  229  to  946  700  yr.  The  evaluation  of  the  periodic  terms  was  carried  out  using  MathCAD  2.0  on  a  Zenith  248PC.  The 
motion  over  5  X  106  yr  shows  that  the  main  terms  combine  to  form  a  quasiperiodic  function  with  a  period  of  about  1. 17  X  106 
yr  and  an  amplitude  of  0*1 5.  However,  over  the  short  term  the  combined  effect  of  these  63  terms  is  an  additional  “precession” 
of  only  1"00X  10  ~  6  yr' 1  which  is  far  too  small  to  be  significant. 

Finally,  cubing  the  third-order  sine  term  in  Eq.  (77)  and  simplifying  using  //Math  a  simple  algebraic  manipulation 
program  on  a  Zenith  248  PC  gives 

—  a  sin  +  Hj )^  =  —  a  sin  e*-^-  |x|^  ^  ^  j  [3  sin(/q/  +  fij)  —  sin(3/qt-|-  3fij)  ] 

.  £  \  (sin  [  (*■.  +  2Kk  )t  +  Hj  +  2fik  ] 

k  Kj  \  Kk  ) 


3_ 

4 


+  sin[(/q  -lKk)t+Hj  -2/xk]  -  2  sin(*/ +  //,  ))( 1  —  8Jk) 

3  toJG}lBmNjNlNm 

+  — XX - - (sm[(/f.  +K,  +Km)t  +  [Xj  +Ht  +//m 

2  lm  KjK,Km 

+  sin[(/r;  -k,  +Km)t  +  Hj  -p,  +  //m]  +sin[(«)  +  k,  -  Km)t 

+  fJ-j  +Vi  -^m]  +  sin[ (Kj  —  k,  -Km)t+fij -fi,  -//m])JJ  , 


(84) 


where  k  is  summed  from  1  to  8,  /  is  summed  from  j  -f  1  to  8,  m  is  summed  from  k  +  1  to  8,  and  Sjk  is  the  Kronecker  delta. 
Equation  (84)  consists  of  980  terms.  Integrating  this  equation  with  respect  to  time  gives 

•  *  1 
a  sin  e*  — 


-  t  ?  “*[<«.+ i 

[(Kj  -2 Kk)t  +  Pj  -  2jik  ]  -  —  cos(Kjt  +  fij) J(1  -8 

Kj  J 

( - l— 

\Kj  +K,+K, 


Kj  ~  2Kk 


3  XX  ajm,G3mNiN,Nm 

2  lm 


KjK,Km 


Xcosf  {Kj  —  K,  +  Km  )t  +  fj.j  -fl,  +flm  \  + 


COS[(A/  +K,  +Km)t+Hj  +fl,  +fim]  + 

1 


1 


Kj  +K,—K, 

+ - - - cos  [  (Kj  -k,  -Km)t  +  Hj -n,  -//„])  ]  +C. 

Kj  Kj  Km  )\) 


Kj  ~K,+K„ 
COS  [( Kj  +K,  -Km)t  +  Hj  +Hl 


(85) 


-20.00  — I - 1 - 1 - 1 - 1 - - 1 - 1 - 1 - 1 - 

-5000000.00  -2500000.00  0.00 


Yean  from  present 

Fig.  6.  Very  long-period  nutation  in  longitude.  The  first-order  change 
in  the  longitude  of  Mars’  pole  with  respect  to  its  orbital  plane  resulting 
from  the  change  in  Mars’  orbital  plane  with  time. 


Evaluating  Eq.  (85),  also  using  MathCAD,  shows  that  the 
amplitudes  of  the  individual  terms  are  as  large  as 
0'.'002  097  +  0"000  006  and  periods  ranging  from  23  807  yr 
to  1.6  X 107  yr.  Evaluation  over  5000  yr  shows  a  change  in 
orientation  of  about  0"003  with  no  discernible  period.  There 
does,  however,  appear  to  be  a  “short-period”  term  with  a 
period  of  about  130  000  yr  and  an  amplitude  of  about 
0"0003.  The  short-period  (  <  10  000  yr)  change  in  orienta¬ 
tion  with  time  is  negligible.  Therefore,  the  third-order  ap¬ 
proximation  of  the  sine  function  does  not  make  a  significant 
addition  to  the  precession  of  Mars. 

APPENDIX  E:  PLANETARY  TORQUES  ON  MARS 

In  the  zeroth-order  estimate  for  the  effect  of  the  planets 
upon  Mars  it  will  be  assumed  that  all  the  planets  in  the  solar 
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Table  6.  Coefficients  for  the  B  J/2  term  for 
the  precession  caused  by  the  planets  on 

Mars. ^ 

Table  7.  The  B  5n/2  coefficients  for  the  nutation 
caused  by  the  Earth  and  Jupiter  on  Mars. 

Coefficient 

n 

Earth 

Jupiter 

Mercury 

2.992 

0 

38.006 

3.418 

Venus 

8.411 

1 

36.203 

2.163 

Earth 

38.006 

2 

32.188 

1.051 

Jupiter 

3.418 

3 

27.254 

0.449 

Saturn 

2.344 

4 

22.265 

0.178 

Uranus 

2.080 

5 

17.701 

0.067 

Neptune 

2.032 

6 

13.775 

0.024 

7 

10.538 

0.009 

8 

7.949 

0.003 

two  assumptions  reduce  Eqs.  (21)  and  (22)  to 


—  P—Q  sin  2A  sin  e, 

2/ 

—  A  sin  e  =  Q  SSL  sin  2e  cos2  A.  (86) 

2  r 

Expressions  for  the  “ecliptic”  longitude  and  distance  ras  a 
function  of  time  need  to  be  derived  before  Eqs.  (86)  can  be 
solved.  Basic  vector  addition  represents  the  position  of  a 
planet  with  respect  to  Mars  by 

r  sin  A  =  rp  sin  Lp  —  rM  sin  St, 

r  cos  A  —  rp  cos  Lp  —  rM  cos  St, 

(87) 

where  rp  is  the  distance  of  the  planet  from  the  Sun,  rM  is  the 
distance  of  Mars  from  the  Sun,  Lp  is  the  heliocentric  longi¬ 
tude  of  the  planet,  and  St  is  the  heliocentric  longitude  of 
Mars.  Under  the  assumption  of  perturbation  free,  circular 
orbits,  rp  and  rM  are  constants  and 

LP  =  npt  +  Lp o, 

St  =  nMt+St  0,  (88) 

where  np  and  nM  are  the  mean  motion  of  the  planet  and 
Mars,  respectively,  and  L ^  and  St  0  are  the  heliocentric  lon¬ 
gitudes  of  the  planet  and  Mars  at  some  epoch  t0.  For  ease  of 
notation,  the  heliocentric  longitudes  of  epoch,  Lp()  and  St  0, 
will  be  dropped  until  they  are  needed  for  evaluation  of  the 
precession  and  nutation. 

Using  Eqs.  ( 87 )  for  sin  A  and  cos  A  in  Eqs.  ( 86 )  gives 


1  ft2  t 

—  p  =  Q  —  sin  e  --  —  sin  In  t  +  — —  cos  lnMt 
r  r  \  2  2 

-rprM  sin(np  +nM))  (89) 

and 

—  A  sin  e  —  Q  SSL.  sin  2e  -1  f—  +  -SL  +  —  cos  In  t 

2t  t  \  2  2  2 


+  — cos  2 nMt  +  rprM  cos (np  +  nM)t 


+  rprMcos(np  -  nM)t 


)■ 


(90) 


The  determination  of  the  value  of  r  as  a  function  of  time 
begins  with  the  cosine  law 

r2  =  tp+r1M-2rprMcos(np-nM)t.  (91) 

In  both  ( 89)  and  (90)  r  appears  as  some  power  r  ~ ',  where  .s 
is  a  positive  integer.  Using  Eq.  (91 )  r~ s  is  rewritten 


r  s=[i2p+tM-2rprMcos(np-nM)t\  5/2 

1  +  (— )  —2  — cos (n  -nM)t 
\rp  )  rp 


(92) 


or 


■COS  (np  —  nM)t 


-s/2 


(93) 


Both  of  Eqs.  (92)  and  (93)  have  the  general  form 
(1  +  x)~s/2,  where  x  =  b 2  —  2b  cos(np  —nM)t  and  b  is 
either  the  ratio  rp/rM  or  rM/rp.  The  solution  to  these  equa¬ 
tions,  provided  the  smaller  of  the  ratios  rp/rM  or  rM/rp  is 
always  used,  is 


+  '2,Bi2cos[k(np-nM)t] 

k 


(94) 


where  ra  is  the  smaller  of  rp  or  rM  and  B  %2  are  the  Laplace 
coefficients.  The  Laplace  coefficients  frequently  arise  in  clas¬ 
sical  celestial  mechanics  multibody  problems  such  as  the 
perturbation  of  the  orbit  of  one  planet  by  others.  The  method 
of  solving  for  the  value  of  any  Laplace  coefficients  can  be 
found  in  several  places  including  Smart  ( 1953 ) .  The  values 
for  B  g/2  for  each  of  the  planets  torque  on  Mars  are  given  in 
Table  6  and  the  values  of  B5/2  through  B{/2  for  the  Earth 
and  Jupiter  with  Mars  are  given  in  Table  7. 

Substituting  Eq.  (94)  into  the  equation  for  the  motion  in 
longitude  of  the  pole  of  Mars,  Eq.  (90),  gives 
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—  X  sin  e  =  sin  B o/2^~(  1  +  a2)  +  cos  2 nMt  -f  -i-  a2  cos  2npt 

+ -2-cos  («p  +nM)t  +  -2-cos  (np  + V5f(-|{  1  +  a2)cos  At(«p  -nM)t 

2  2  /  k  \  2 

1  1  ^ 

+  —  cos  [Acnp  -  (k-2)nM]t  +  —  cos  [At«p  -  (Ac  +  l)nM  ]t  +  -L  {cos[  (k  -  2)np  -  knM]t 

+  cos[  (Ac  +  2)np  -AcnM]f}  +  ^a{cos[(Ac  +  1  )np  -  ( k -  1  )nM]t 

+  cos[(Ac-  l)«p  -  (At-f  \)nM]t  +  cos[(Ac  +  \){np  -  nM)  ]t  +  cos[  (Ac  +  1 )  («p  .  (95) 


This  equation  is  integrated  to  yield 


—  A  sin  e  =  sin  2e 


(7'“ 


11a2 
—  ( 1  +  a2)t  +  — - sin  2 nMt  +  — —  sin  2 npt 


8  nM  "  8np 


2  U»+«. 


1  sin(«p  +  nM)t  + - — 


-  sin  (n 


—  «Af)^  j  J 


+  If^5/2 

k  2 


(1  +  a2) 


k(n.  -  nM ) 


■sin  Ac(«p  —  nM)t  + 


2 [knp  -  (k-2)nM] 


1 


8[Acwp  -  (k  +  2)nM] 
2  f  1 


+  a 

+ 


2  l  [ (A:  —  2)np  -AcnM] 
1 


sin [ Ac«p  -  (k  +  2)nM]t 


sin[(Ac-2)np  -  knM]t  + - 


[(Ac+  l)np  -  (Ac-  1  )nM] 
1 


[  (A:  +  2)np  —  knM] 
cos[ (A:  +  l)np  -(k-\)nM]t 


sin [ knp  -  (k-2)nM]t 


sin [ ( Ar  -+-  2)np  -  knM]t 


[(At-  l)np  -  (At 4-  1  )nM] 
1 


sin [ (At  —  l)«p  -  (At+  1  )nM~\t 


[  (At  +  l)(«p  —  nM)] 

_ 1 

[(At+  l)(«p  -«M)] 


sin[ (At -f  l)(«p  -  nM)]t 

[(At+  l)(«p -nM)]t J  J  . 


-cos 


(96) 


Like  any  equation  involving  numerous  sine  functions  Eq.  (96)  looks  daunting  but  it  is  actually  quite  simple. 
Similarly,  substituting  the  equation  for  the  expression  r~s/ 2  into  the  equation  for  the  motion  in  latitude  gives 


 Q» 


in  f  j/i.i,![a- 


sin  2 npt  +  sin  2 nMt  —  2 a  sin(«p  —  nM)t  ] 


+  'ZBk2(<^r{sin[(k  +  2)np  -  Ac«M]H- sin[  (Ac  -  2)np  -  knM]t} 

k  V  2 

+  i  {sin [ knp  -  (Ac  -  2)nM]t  +  sin[Atnp  -  (k  +  2)nM]t} 

-a{sin[(Ac  +  l)«p  -  (At  -  \)nM ]t  4-  sin[  (Ac  -  l)np  -  {k+  l)«M]f}jj  . 


(97) 


The  integral  of  Eq.  (97)  is 


P  ="TT_sin  £\Bo/2(-~ —  cos  2«pt  +  — ^ —  cos  lnMt - — — Cos(np  -  «M)t) 

2  ra  l  \2  np  2  nM  np-nM  ) 


(At-f  2)np  -knA 


-cos[(Ac  +  2)np  -knM]t  + 


( Ac  —  2 )  «p  -  Ac«a 


-cos  [(At  —  2)«p  —  knM]t 
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+  t(t - 71 — ^ - cos [knp  -  ( k-2)nM]t  +  -  1  - cos [knp  -  (k  +  2)nM]t\ 

2  \  knp  -  (k  -  l)nM  knp-(k+2)nM  J 

-i 


+ 


(&+  1  )np  —  (k—  1  )nM 
1 


cos[(fc+  \)np  -  (k-  1  )nM]t 


(k  —  l)n  —  (k  +  1 


y— -cos [(k-  \)np  -  (k+  \  )nM]t  jjj 


(98) 
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